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ABSTRACT 

Several UV and near-infrared color selection methods have identified galaxies at z = 1-3. Since each 
method suffers from selection biases, we have applied three leading techniques (Lyman break, BX/BM, 
and BzK selection) simultaneously in the Subaru Deep Field. This field has reliable (Az/(l+z) = 0.02- 
0.09) photometric redshifts for ~53,000 galaxies from 20 bands (1500 A-2.2 /jm). The BzK, LBG, 
and BX/BM samples suffer contamination from z < 1 interlopers of 6%, 8%, and 20%, respectively. 
Around the redshifts where it is most sensitive (z ~ 1.9 for star-forming BzK, z ~ 1.8 for z ~ 2 
LBGs, z ~ 1.6 for BM, and z ~ 2.3 for BX), each technique finds 60%-80% of the census of the 
three methods. In addition, each of the color techniques shares 75%-96% of its galaxies with another 
method, which is consistent with previous studies that adopt identical criteria on magnitudes and 
colors. Combining the three samples gives a comprehensive census that includes ~90% of z p hot = 13 
galaxies, using standard magnitude limits similar to previous studies. In fact, we find that among 
z = 1-2.5 galaxies in the color selection census, 81%-90% of them can be selected by just combining 
the BzK selection with one of the UV techniques (z ~ 2 LBG or BX and BM). The average galaxy 
stellar mass, reddening, and star formation rates (SFRs) all decrease systematically from the sBzK 
population to the LBGs. and to the BX/BMs. The combined color selections yield a total cosmic SFR 
density of 0.18±0.03 M Q yr" 1 Mpc~ 3 for K AB < 24. We find that 65% of the star formation is in 
galaxies with E(B — V) > 0.25 mag, even though they are only one-fourth of the census by number. 

Subject headings: galaxies: distances and redshifts - galaxies: evolution - galaxies: high-redshift 
galaxies: photometry - infrared: galaxies - ultraviolet: galaxies 



1. INTRODUCTION 

Several color techniques have succeeded in identifying 
large samples (thousands) of galaxies in various windows 
of high rcdshift. They work by using deep wide- field 
imaging in only a few broad-band filters. In the sim- 
plest cases, only two colors are needed to isolate a spec- 
tral break. For imaging surveys limited to optical ob- 
servations, the Balmer/4000 A break is measurable up 
to z ~ 1, and the Lyman continuum break is detectable 
starting at redshifts of 3 and higher (e.g., Steidel et al. 
1999; Bouwens et al. 2006; Yoshida et al. 2006). But 
at intermediate redshifts, CCDs are only sensitive to 
the spectral region between these two strong features. 
The resulting inability to identify galaxies at z w 1-3, 
has lead to this range being called the "redshift desert." 
Finding large samples of galaxies in the redshift desert re- 
quires detections of the Balmer break with near-infrared 
photometry, or the Lyman break with ultraviolet imag- 
ing. 

The extension of the Lyman break technique to z < 2.6 
requires deep near-ultraviolet (NUV) imaging, which has 
only been recently available with the GALEX satellite 

9 Current Address: Space Telescope Science Institute, Balti- 
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(Ly et al. 2009, hereafter L09) and Hubble WFC3/UVIS 
(Hathi et al. 2010). Prior to this, different techniques 
were developed, to select "BX" and "BM" galaxies (Adel- 
berger et al. 2004; Steidel et al. 2004), BzK galaxies 
(Daddi et al. 2004), and "distant red galaxies" (DRGs; 
Franx et al. 2003; van Dokkum et al. 2004). However, 
each technique suffers from its own selection biases. For 
example, UV selections (e.g., LBG, BX) tend to iden- 
tify young star-forming galaxies with low dust extinc- 
tion, while IR techniques (e.g., BzK, DRG) select more 
massive, dusty galaxies. 

In this paper, we apply several color selection tech- 
niques to our panchromatic photometry of the Subaru 
Deep Field (SDF), to identify «21,000 galaxies at z = 1- 
3 for a census of optical and near-infrared (NIR) selected 
star-forming galaxies. Throughout this manuscript, the 
term "census" will be repeatedly use to refer to the 
union of the different color selection methods that iden- 
tify LBGs, BX/BMs, and sBzK galaxies down to specific 
magnitude depths that are typical of past and current 
z ~ 2 census studies. Our multi-technique survey, com- 
pared to previous work (Rcddy et al. 2005; Quadri et al. 
2007; Grazian et al. 2007, hereafter R05, Q07, and G07, 
respectively) , has significantly more galaxies because our 
deep imaging covers 2-10 times more area. Since we 
also have a wide range of photometry spanning many 
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filters, we can then directly compare these simple color 
techniques to a sample derived from photometric redshift 
(Baum 1962, hereafter photo- z or z p hot)- 

The observations spanning wavelengths from 1500 A to 
4.5 /urn and their reductions are described in Section 2. 
Section 3 describes the merging of the multi-wavelength 
data together, which accounts for differences in the spa- 
tial resolution across different wave bands. Section 4 dis- 
cusses the determination of photo- z and the modeling of 
the spectral energy distribution (SED). In Section 5, the 
different color selection techniques used and individual 
sample properties are presented. Section 6 discusses the 
selection effects of one technique against another, and 
compares them with our z p hot-selected sample. We also 
compare our results with previous surveys. In Section 7, 
we present a Monte Carlo simulation that aims to gener- 
ate a mock census of z k 1-3 galaxies to support many 
of our observed relations and results. Conclusions from 
our z — 1-3 census survey are summarized in Section 8. 

Throughout this paper, a flat cosmology with Hm, 
^70] = [0-7, 0.3, 1.0] is adopted. Magnitudes are reported 
on the AB system (Oke 1974), and unless otherwise in- 
dicated, limiting magnitudes are 5cr, corrected for the 
flux falling outside of the apertures assuming an unre- 
solved source (i.e., total limiting magnitudes). Diameter 
apertures sizes are denoted with a u <f>." 

2. SUBARU DEEP FIELD MULTI-WAVELENGTH 
OBSERVATIONS 

A summary of the depth, spatial resolution, and sam- 
ple sizes for the imaging data is provided in Table 1 , while 
the spatial coverage of the SDF at different wavebands is 
shown in Figure 1. Although the reported depths illus- 
trate the sensitivity of the data, we emphasize that all 
of the sources actually analyzed later in this paper are 
detected well above these 5a limits in almost all of the 
wavebands. The only exceptions are the FUV, NUV, J, 
H, and K bands. Since most of these data are described 
in previous papers, only brief summaries are presented. 

2.1. Subaru/ Suprime- Cam Data 

The BVRci'z' observations 3 are the deepest data ob- 
tained to date from Suprimc-Cam (Miyazaki et al. 2002), 
and are described in Kashikawa et al. (2004). The total 
limiting magnitudes were redetermined with more rig- 
orous masking of pixels affected by object flux and arc 
« 0.1-0.2 mag deeper than those quoted in Ly et al. 
(2007). They range from 26.1 to 27.9 mag. 

In addition to these data, Suprime-Cam imaging in 
four intermediate-band filters, IA598, IA679, z&, and z r 
are also available. Studies (e.g., Meisenheimer & Wolf 
2002; Ilbert et al. 2009; van Dokkum et al. 2009) us- 
ing intermediate-band filters have improved photo-z es- 
timates, as they are able to identify the location of spec- 
tral breaks more accurately. The Zf, and z r observations 
were obtained in 2003 and 2004, and are described in 
Shimasaku et al. (2005). The net integration times are 
85 and 600 minutes for Zf, and z r , which give total limit- 
ing magnitudes of 25.6 and 24.7 mag, respectively. The 
IA598 and IA679 observations were obtained in 2007 and 

3 These data are publicly available at 
http : / /soaps . nao . ac . jp/SDF/vl/index.html. 
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Fig. 1. — Sky coverage of SDF data described in Section 2. 
Masked regions for the Suprime-Cam data are shown with the 
solid black rectangles and circles. Excluding these regions, the 
area coverage is 870.4 arcmin 2 . The shaded regions correspond to 
lower i^-band sensitivity areas (718.7 arcmin 2 has deep K cover- 
age). GALEX, NEWFIRM, and MOSAIC data have coverage over 
the entire SDF. Spitzer 1 ks IRAC coverage is illustrated with the 
dashed line (blue; [3.6] and [5.8]) and dot-dashed line (green; [4.5] 
and [8.0]). The angular size for 2 hy Mpc is shown by the thick 
bars for z = 0.5 and z = 1-3. (A color version of this figure is 
available in the online journal.) 

are described in Nagao et al. (2008). Their depths are 
26.0 and 26.7 mag, respectively. 

Finally, we include imaging in five narrow-band (NB) 
filters (NB704, NB711, NB816, NB921, and NB973). 
These data are included to specifically help improve the 
Zphot's for emission-line galaxies, but we find redshift ac- 
curacy improvements for galaxies at z < 1.6. Their 5cr 
total flux depths vary between 24.9 and 26.2 mag. The 
narrow-band data reduction is further discussed in Ouchi 
ct al. (2003), Shimasaku et al. (2003), Shimasaku et al. 
(2004), and Kashikawa et al. (2004). 

2.2. UKIRT/WFCAM Data 

K-h&ml data for this survey were acquired with Wide- 
Field Camera (WFCAM; Casali et al. 2007) on the 
United Kingdom Infrared Telescope (UKIRT) on 2005 
April 14-15 and 2007 March 5-6. Because the field of 
view (FoV) of WFCAM is split into four 13f6xl3.'6 im- 
ages separated by 12. '8, four pointings were required to 
cover the whole SDF optical image. The total exposure 
time for each pointing is 294-300 minutes, except for one 
pointing of 15 minutes. We ignore the less sensitive re- 
gions, so 80% of the SDF optical coverage has uniform K- 
band data down to the deepest sensitivity (see Figure 1). 
Previous results using a subset of this data were pre- 
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TABLE 1 

Summary of SDF Multi-wavelength Data 
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303 


0.91 


26.276 
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26.023 
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1 Filter's center and FWHM are in units of A. 

^ Limiting magnitudes arc defined as 5(7 with a 3" cj> (^j J i an d H), 6'.'8 (ft (FUV and NUV), or l" (ft (remaining wave 
bands) aperture. IRAC limits are provided with a A" and 2" (j) for original and HiRcs data products. 

^ Corrections of aperture flux to total source flux assuming a point source. These are determined within the same aperture 
as mj^. 

4 The 5<r aperture-corrected limiting magnitudes. 

^ First values exclude sources in the if -band low sensitivity regions. Those in parentheses include sources in these regions. 

0.06 mag) in the predicted U 



sented in Hayashi et al. (2007). The data were reduced 
following standard NIR reduction procedures in IRAF. 
Photometric and astrometric calibrations were carried 
out against Two Micron All Sky Survey (2MASS; Skrut- 
skie et al. 2006) stars. The final if-band data reach a 
total ha depth of 23.5 mag. 

2.3. Mayall/MOSAIC Data 

[/-band data were obtained from the KPNO Mayall 4 
meter telescope with the Mosaic-1 Imager (FoV of 36'; 
Muller ct al. 1998) on 2007 Apr 18 and 19. Observing 
conditions were dark with minimal cloud coverage. A se- 
ries of 25 minute exposures was obtained to accumulate 
47 ks for a 5<r total depth of 26.1 mag. A standard dither- 
ing pattern was followed to provide more uniform imag- 
ing across the «10" CCD gaps. The MSCRED IRAF 
(version 2.12.2) package was used to reduce the data and 
produced final mosaic images with a pixel scale of 0'.'258 
and an average- weighted seeing of wl''5 FWHM. The re- 
duction steps followed closely the procedures outlined for 
the reduction of the MOSAIC data for the NOAO Deep 
Wide-Field Survey. 

While several Landolt (1992) standard star fields were 
observed, Landolt's U filter differs significantly from the 
MOSAIC/?/, which makes it difficult to calibrate the 
photometry. One hundred and two Sloan Digital Sky 
Survey (SDSS) stars distributed uniformly across the 
eight Mosaic-1 CCDs with v! < 21 mag were used. This 
approach requires a transformation between SDSS vl and 
MOSAIC/f/, which is obtained by convolving the spectra 
of 175 Gunn-Stryker stars with the total system through- 
put at these wavebands. These stars span B — V = —0.2- 
0.8 mag, and show that the B — V color term is smaller 



than the scatter (la 
magnitudes. 



2.4. GALEX/NUV and FUV Data 

The GALEX space telescope (Martin et al. 2005) pro- 
vides simultaneous FUV and NUV imaging with a pixel 
scale of 1'.'5 and an FoV of 1?2, thus covering the entire 
SDF. A total of 161225 s 4 (80851 s) was obtained in 
the NUV (FUV) band (Malkan 2004). The data were 
processed through the GALEX reduction pipeline and 
stacked accordingly. The FWHM of unresolved images 
is 4'.'5 in the FUV and 5'.'0 in the NUV. The 5cr total 
limiting magnitudes are both «25.9 mag. 

2.5. Mayall/NEWFIRM Data 

J- and ii-band observations were obtained from the 
NOAO Extremely Wide-field Infrared Imager (NEW- 
FIRM; Probst et al. 2004, 2008) on 2008 May 10-12. 
Exposures of 60 (J) and 20 or 30 (H) s were taken with 
a simple five-point dithering pattern to image across the 
detector gaps (« 35" in size) and to avoid latent images 
of bright sources on the same source. A total of 27.5 
ks and 5.5ks were acquired in J and H, respectively. 
While cloud coverage was minimal, the conditions were 
non-photometric (e.g., the effective J exposure time was 
approximately 9ks). These NEWFIRM data required a 
first-pass reduction to create a deep object mask. This 
mask is then used to create flat-fields and sky images to 
enable proper subtraction of the sky background. Then 
the sky-subtracted images were used for the second-pass 

4 The results of L09 used a data set that is shallower by 20 ks. 
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processing to create the final mosaics. The J (H) limit- 
ing magnitude was determined to be 23.4 (22.8) mag in 
a 3" <fi aperture, and the average seeing is ~l'/3 (1'.'2). 

2.6. Spitzer/IRAC Data 

Imaging in the IRAC bands was conducted in 2006 
February 11 and 14 (Prop#20229; Malkan et al. 2005). 
The observations consisted of nearly uniform coverage 
in all four bands of 1 ks. The coverage is 26' x 25'.5 
for [3.6] and [5.8] or [4.5] and [8.0], and the overlap- 
ping region for all four bands is approximately 19' x 25f5. 
Version S13.2 data products were used. Necessary pre- 
processing steps were taken to remove detector arti- 
facts (e.g., mux-blccding, column pull-downs) with con- 
tributed software provided by the IRAC community, and 
the data reduction followed standard procedures using 
the MOsaicker and Point source Extractor (MOPEX) 
package. The astromctric solutions for the mosaics were 
examined against 2MASS stars and are accurate to with 
071. 

Apertures were placed randomly on the images to de- 
termine the limiting depth. The depths are summarized 
in Table 1. For the remainder of this paper, we only 
use the [3.6] and [4.5] data, since they sample rest frame 
~ 1-2 fim for z = 13, which arc unaffected by polycyclic 
aromatic hydrocarbon (PAH) features. 

The greatest limitation of the IRAC data is the spa- 
tial resolution: measuring unresolved galaxies generally 
requires at least a 4" <f> aperture. Recently, Velusamy 
et al. (2008) developed a high-resolution image decon- 
volution code (HiRes) to handle Spitzer data. The code 
has been found to significantly improve the quality of the 
images, allowing for smaller photometric apertures and a 
reduction of source confusion. HiRes was used on these 
IRAC data, and the results show improvements in the 
aperture fluxes by a factor of 2-4 using the same aper- 
ture sizes. For example, a 2" (f> aperture was able to en- 
compass 95% and 81% of the total light in [3.6] and [4.5], 
respectively. We examined sources in both images and 
found that photometric fluxes are conserved, indicating 
that these HiRes images can be used to attain reliable 
flux measurements in these bands. This improvement al- 
lowed for the crucial detection of lower mass galaxies at 
z ~ 2 at rest 1 /xm. 

3. THE SDF MULTI-BAND GALAXY SAMPLE 
3.1. Synthesizing Multi-band Data 

Methodology. A key requirement of this survey is to 
merge the full SDF multi-wavelength data. However, the 
spatial resolution varies by a factor of a few for the differ- 
ent instruments. Simply degrading the best data (~1") 
to match the poorest resolution (4"-5" for GALEX data) 
will hamper the identification of faint galaxies, most of 
which are at high redshift. 

A solution to this problem is a hybrid approach of 
(1) slightly degrading high-resolution images to match a 
"baseline" spatial resolution that is not significantly dif- 
ferent and (2) using aperture corrections. This method 
will not greatly compromise sensitivity, and the aperture 
corrections are reliable, since this work mostly probes 
compact galaxies (they are faint and at high-z). 

We begin by combining multiple broad- and 
intcrmcdiate-band images into a single ultra-deep frame. 



This approach benefits from the inclusion of weak sources 
that are marginally detected in individual wave bands, 
but are highly significant in the merged image. This 
ultra-deep image serves as a reference image when run- 
ning SExtractor (Bertin & Arnouts 1996) in "dual im- 
age" mode. It yields a catalog of sources in the SDF, 
and later serves as the "master object list" for the photo- 
z catalogs and the color-selected samples. The similari- 
ties of the spatial resolution for BVRci' 'z'K , IA598, and 
IA679 allow us to degrade these images to a l'.'l (here- 
after "psf-matched" ). The inclusion of K implies that 
this ultra-deep mosaic is simultaneously sensitive to op- 
tically selected and NIR-selected galaxies. 

To stack these images, we use the x 2 -weighting scheme 
of Szalay et al. (1999). This method uses the x 2 of each 
pixel, imaged in N bands, to estimate the probability 
that it is sampling a source or the sky. First, for image i, 
the mean (/z^) is removed, and then the image is scaled 
by the sky rms fluctuation (<Tj): 

9i(x,y) = , (1) 

where fi(x,y) is the measurement (in count rates, elec- 
trons, or fluxes) in image i, and g%{x, y) can be viewed as 
the signal-to- noise ratio (S/N) for a given pixel. Then, 
we sum the squares of these 8 S/N images to yield the 

8 

ultra-deep image: y = g\ . 

i=l 

We considered including additional wave bands into 
the ultra-deep image, but they degraded the results due 
to their lower sensitivities and/or poorer spatial resolu- 
tions. For example, the U-, J-, and H-band data have 
lower spatial resolution (~ 1'.'5) and are less sensitive 
than the nearby B- and K-hand data. Similarly, includ- 
ing the narrow-band observations does not increase the 
sample sizes. These data are shallower than the broad- 
band SDF observations, so almost all of the narrow- 
band selected sources are identified in the broad-band 
images. While there are some galaxies that are only de- 
tected in the narrow-band observations and thus missed 
by our multi-band x 2 image (e.g., Lya emitting galaxies 
at z > 5), this tiny galaxy population is less than 0.05% 
of the full catalog, and will not affect the z = 1-3 galaxy 
census that we later derive. For these Lya emitters, it 
is better to generate a separate sample using only the 
narrow-band observations, as done in Kashikawa et al. 
(2006) and Shimasaku et al. (2006), for example. 

Source catalogs. Wc run SExtractor on the square root 
of the ultra-deep \ 2 image to obtain a total of 243,964 
sources with a minimum of five connected pixels above 
a x — 3.5 threshold, of which 211,594 sources are in 
the unmasked regions shown in Figure 1. This sample 
is further reduced to 174,837 sources when the poorly 
sampled if-band regions are excluded. This is the largest 
photometric catalog for the SDF. 

One concern with the multi-band catalog is that it may 
systemically miss some galaxies, particularly the bluest 
and reddest ones, because they are undetected in the 
longer or shorter bandpasscs, respectively. We inves- 
tigate this possibility by generating BVRc and i' z'K 
stacked mosaics using the x 2 weighting approach. Wc 
run SExtractor on these stacked images, and compare 
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these catalogs against the ultra-deep eight-band catalog. 
First, we find that the overlap is «90% or larger at any 
given magnitude in any of the individual bands. Second, 
since galaxies with extreme colors could be systemati- 
cally missed, we also examine the overlap of the bluest 
(reddest) quartile of BVRc (i'z'K) catalog against our 
eight-band catalog. The blue sample consists of —20,000 
galaxies with B — V < 0.2 mag while the red sample 
contains —6000 galaxies with z' — K > 0.4 mag. For- 
tunately, the overwhelming majority of these galaxies 
with extreme colors were already included in the origi- 
nal selection using the eight-band ultra-deep image, 93% 
and 89% of the blue and red samples, respectively. We 
find that these fractions hold regardless of magnitudes 
or colors. These tests and comparisons indicate that the 
multi-band synthesis approach does not inherently miss 
any significant population of galaxies. We will discuss 
this issue further in Section 7. Szalay et al. (1999) also 
find that the approach is optimal when generating source 
catalogs for the Hubble Deep Field. 

Measured fluxes. Generally, fluxes are measured in 
fixed circular apertures with diameter sizes comparable 
to twice the FWHM of resolution, and are then corrected 
for aperture losses, which are listed in Table 1. Thus, 
for images with < 1", we measure photometric fluxes 
in 2" (f> aperture with SExtractor. For the U, J, and H 
bands, we use 3" cf> aperture measurements. The GALEX 
and IRAC/HiRcs bands use a 6'.'8 and 2" <j> apertures, 
respectively. For the GALEX and IRAC/HiRcs data, 
the images were not matched to the Suprime-Cam pixel 
scale. Instead, we use DAOPHOT to measure the fluxes 
at known positions, which follows the identical procedure 
discussed in L09. One advantage of DAOPHOT is that 
additional noise is not introduced when these images are 
regridded to the reference SDF Suprime-Cam image with 
a pixel scale of 0'.'2. 

GALEX source confusion. One concern for GALEX 
FUV and NUV measurements is source confusion — the 
overestimation of the flux for faint sources that hap- 
pen to fall in the wings of brighter nearby sources. To 
correct for this, we take a point-spread function (PSF)- 
fitting approach, which supersedes the baseline approach 
of fixed aperture measurements. This consists of two 
stages: (1) fitting and removing relatively bright sources 
with an empirical PSF constructed from a dozen unre- 
solved bright sources, and (2) then for fainter sources, 
aperture photometry can be done on the image with the 
relatively bright sources removed. The GALEX fluxes for 
the bright and faint sources are then combined together 
for a single catalog. 

The PSF fitting is done using a few DAOPHOT tools. 
DAOPHOT assumes that sources are unresolved relative 
to the PSF. Since 97% (51%) of the objects in our multi- 
band source catalogs are 3" (l'.'l) FWHM or smaller, the 
large spatial resolution of 4"-5" for GALEX implies that 
modeling and removing sources with a point-like PSF is 
a valid approach for almost all sources. 

Prior inputs to fitting the GALEX observations are 
based on detections in the U band. A total of 43,839 
(31,084) sources were considered for NUV (FUV) PSF 
fitting down to 26.75 mag (for NUV) and 26.25 mag 
(for FUV). which is much deeper than any previous 
photometric studies with GALEX. Only 18,036 (14,341) 
sources were actually PSF-fittcd, since the remaining 



sources were either too faint or too close to another 
brighter object. 

While PSF fitting is done for bright sources, we still 
determine the total fluxes using the fixed aperture flux 
measurements and apply aperture corrections. This en- 
sures that all fluxes are obtained in a consistent man- 
ner. We note that we compare "total" flux measurements 
with those obtained from PSF-fitting for bright sources, 
and find them to be consistent within the uncertainties. 
This is expected, since the aperture corrections are de- 
rived from the empirical PSF, which is used to fit and 
remove the flux from the wings of sources. One form 
of "quality assurance" is to compare the "noise" in the 
original GALEX image to that of the "cleaned" image, 
after PSF-fitting subtraction. Visual inspection shows 
that the source-subtracted image looks very much like 
pure white sky noise. We quantified this by measuring 
the la fluctuations in the same area with same clipping 
algorithm, and find that the rms is reduced by 50% after 
the sources are PSF-fitted and removed. In fact, the rms 
of the cleaned NUV image, 19 photons per pixel, is only 
slightly larger than the predicted rms based on Poisson 
count statistics (~ 16). Similarly, the cleaned FUV im- 
age has an average sky background of 13 photons per 
pixel. The square root of this is 3.6 photons, which is 
very close to our observed rms fluctuations of 3.9 pho- 
tons. Thus the sky noise in our cleaned picture is nearly 
as small as the limit set by the photon statistics of the 
sky background. 

3.2. The Identification of Stars 

To identify foreground stars we use two techniques. 
First, stars can be easily identified in the NIR, since it 
measures the Raylcigh- Jeans tail, and forms a stellar lo- 
cus in color space distinct from galaxies. Studies that 
have selected BzK galaxies also distinguished stars by 
their z' — K colors (see e.g., Q07). For sources detected 
in the K band, we classify 1603 stars with z' — K < 
0.455(B - z') - 0.773 and z 1 - K < 0.235(5 - z') - 0.279. 

For sources that are undetected in K (or lack K data) , 
we use the technique (called the A method) that is de- 
scribed in L09. First, objects were considered potential 
stellar candidates based on how similar they are to the 
PSF. They were assigned a number between and 10 
(10 being most PSF-like). However, unresolved galaxies 
might also be very compact. Therefore, to distinguish 
these galaxies from stars, we then calculate their devia- 
tion from the stellar locus in the B — V, V — Rc, and 
Rc — z' colors. This deviation (A) from the locus in 
a multi-color space is certainly affected by photometric 
scatter for faint sources, so we limit stellar identification 
with the A method to V = 25.0. A values for -67,000 
sources are illustrated in Figure 2, and show a strong 
peak at A m (the secondary peak is for low-z galax- 
ies). 

One of the best tests of this method is to compare with 
the sample identified from B — z 1 and z' — K colors. The 
distribution of A for the 1603 JC-band selected stars is 
plotted in Figure 2 and supports the A method. Among 
the if-band stellar sources, 82% (1315/1603) of them 
have point-like rank values of 4 and higher. Thus, we 
classify undetected -fsT-band sources with rank values of 
at least 4, V < 25.0, and A values within 3er of the 
expected fluctuations in A-magnitudc relation as stellar 
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sources. This yielded 1167 additional sources, for a total 
of 2770 stellar sources. We removed these stellar sources 
from our full catalog and report stellar contamination for 
individual galaxy selection samples below. 

4. PHOTOMETRIC REDSHIFTS AND SED MODELING 
4.1. Photometric Redshifts 

For this survey, the number of spectroscopic redshifts 
is limited (N ~ 1000), so well calibrated z p hot's are used. 
These Zphot's are determined from 20 bands (FUV, NUV, 
UBVR c i'z'z b z r JHK, IA598, IA679, and five narrow- 
band filters) using the Easy and Accurate z p hot from 
Yale (EAZY; Brammer et al. 2008) package. This photo- 
z code has the following advantages: (1) it is easy 
to use and fast, (2) it handles input measurements in 
fluxes rather than magnitudes (e.g., hyperz; Bolzonella 
et al. 2000), which is crucial for non-detections (see be- 
low), and (3) it has been well calibrated against existing 
optical-to-NIR data with spectra for a z /(l + z) « 0.03. 

Throughout the paper, rcdshift distribution is de- 
scribed by N(z), which uses the peak of the probability 
rcdshift distribution. We have compared N(z) for a sub- 
set of galaxies at high rcdshift to the distribution made 
by adding the full rcdshift distribution for each galaxy. 
We find that the two distributions are in agreement, jus- 
tifying the use of the peak. 

Illustrated in Figure 3 is a comparison of photomet- 
ric and spectroscopic redshifts for 495 sources using 
20 bands. The spectroscopy was obtained from Sub- 
aru, MMT, and Keck. The inclusion of the NUV 
improved photo-z accuracy by ~17% compared to 
UBVR c i'z'z b z r JK, IA598, and IA679 alone. H and 
F(7Fmeasurements added a further ~10% improvement. 
Photo-z's derived from 15 bands (excluding narrow- 
bands) are accurate to 2.5% out to z « 1.8, and «1.2% by 
including the narrow-band measurements. Spitzer/IRAC 
measurements were also included, but produced little im- 
provement. In particular, we found more problems at 
lower redshift, and this is likely due to PAH features that 
have not been included in the photo-z templates. We 
find that the catastrophic failure rate is ~10%. A pre- 
liminary examination indicates that one-fourth of these 
failures are due to galaxies with very strong optical emis- 
sion lines. 

Our photometric redshifts of galaxies with emission 
lines identified with narrow-band filters (see Ly et al. 
2007) are also found to be reliable. Among 5029 narrow- 
band excess line emitters between z = 0.24 and z = 1.47, 
4568 have z p hot estimates. This is almost 10 times the 
number of spectra that we used to compare z p hot with 
Zspcc- We estimate the median z p hot and find good to rea- 
sonable agreement prior to including narrow-band mea- 
surements and additional improvement with the narrow- 
band filters included (see Table 2). 

These z p hot's derived without the inclusion of narrow- 
band data confirm (1) the ability of Ly et al. (2007) to 
distinguish Ha, [O III], and [O II] emitters accurately 
using simple two-color selection, and that (2) the z p hot 
derived for SDF galaxies are reliable. 

Another illustration of the accuracy of the z p hot's is 
shown in Figure 4, where half of the upper 68% minus 
the lower 68% confidence ranges in redshift is shown ver- 



TABLE 2 

Photometric Redshifts for Emission-line Galaxies 



Sample 



2 p hot(15-band) z phot (20-band) 



Ho NB816 
Ho NB921 



O III 
O III 
O III 
O III 



NB704 
NB711 
NB816 
NB921 



O II] NB704 
O II] NB711 
O II] NB816 
O II] NB921 



0.243 
0.401 
0.407 
0.423 
0.630 
0.837 
0.891 
0.912 
1.189 
1.467 



0.257±0.015 
0.410±0.008 
0.409±0.007 
0.428±0.004 
0.635±0.009 
0.788±0.010 
0.797±0.088 
0.890±0.013 
1.093±0.272 
1.537±0.033 



0.233±0.003 
0.396±0.007 
0.403±0.000 
0.431±0.000 
0.645±0.046 
0.835±0.000 
0.890±0.000 
0.909±0.000 
1.193±0.002 
1.473±0.000 



TABLE 3 

Summary of Photometric Redshift Uncertainties 









5BK3cr 




5B5cr 


N 


Median ± a 


N 


Mcdian ± a 


0.0- 


-0.5 


3213 


0.020±0.008 


6763 


0.024±0.014 


0.5- 


-1.0 


7029 


0.019±0.009 


16172 


0.020±0.008 


1.0 


-1.5 


4508 


0.019±0.006 


10205 


0.025±0.013 


1.5 


-2.0 


4242 


0.033±0.016 


12827 


0.062±0.038 


2.0- 


-2.5 


1673 


0.073±0.029 


9257 


0.091±0.029 


2.5- 


-3.0 


1253 


0.044±0.023 


5313 


0.050±0.025 


3.0- 


-3.5 


522 


0.050±0.032 


2388 


0.041±0.021 


3.5- 


-4.0 


353 


0.019±0.005 


1333 


0.018±0.004 


4.0- 


-4.5 


45 


0.020±0.005 


374 


0.019±0.006 


4.5- 


-5.0 


2 


0.017±0.006 


30 


0.014±0.002 


NOTE. 


— Reported arc the median values 


and standard 



deviations of the 68% z p hot uncertainty: Az68 

Z_6s)/2/(l + Zphot)- 



(z+6 



sus photo-z. We refer to this as Az 68 e 



(z + 68— Z-6s)/2 
1 + Zpho 



and use it to indicate z p hot accuracy. The median and 
standard deviation of Azes are reported in Table 3 and 
are shown in Figure 4. As expected, z sa 1.8-2.8 suffers 
from larger z p i 10 t errors due to the absence of a Lyman or 
Balmer break within the CCD filters, lower sensitivity in 
the NIR, and in most cases, non-detection in the rest-UV 
wave bands. 

We did a test excluding if -band measurements (our 
deepest NIR data) for a set of 24,000 if-band galax- 
ies detected at 3<7 in at least five bands. We find that 
the z p hot errors increased by 20% (mostly at z > 1), 
that the number of outliers (see below) doubled, and 
the greatest (Az pno t ~ 0.4) systematic offsets occurred 
at Zphot = 1.3-2.2. However, for the majority («75%) 
of sources, the photo-z's without K are good to within 
Az = ±0.05, particularly at z p hot < 1-2. Therefore, 
we generate two photo-z samples consisting of bright 
sources with if-band data and those that lack it. The 
first sample requires a minimum of five-band detection 
at the 3cr limit (hereafter 5BK3cr) with one of the bands 
being K. The second sample consists of at least five- 
band detection (without any K restrictions) at the 5<7 
level (hereafter 5B5ct). Note that only the measurements 
from broad and intermediate bands (i.e., narrow bands 
are excluded) are used to construct these photo-z sam- 
ples. These samples contain 22853 and 64691 galaxies 
(after stellar removal). Combining the non-stellar 5B5<7 
and 5BK3er samples yield 65117 galaxies. Since 15-band 
detection is not available for all sources (e.g., some are 
undetected in the NUV), we summarize the mixture of 
detections that we have across 15 bands in Figure 4.1. 
Detections in at least 10 wave bands represent 99% of 
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Fig. 2. — Deviation (A) in magnitudes from the stellar locus in a three-color space. (Left) Distribution of A is shown for ~67,000 sources 
with point-like rank of 4—10 in the top panel and against V in the bottom panel. The solid lines correspond to a minimum A = ±0.1 mag 
and the ±3<r fluctuations expected from a combination of the B, V, i?c, and z' depth. (Right) The distribution of A for K-band selected 
stellar sources based on their B — z' and z' — K colors. This shows that the majority of stellar sources have A rs 0. 
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Fig. 3. — Comparison of photometric and spectroscopic redshifts. A total of 495 spectroscopic sources have z p hot estimates. The z p hot 
are derived from 20-band photometry and show accuracy of (z p hot — z spoc )/(l + z Bpcc ) as 0.012 with less than 1% systematic offset. We 
find similar results using a smaller sample that has > 5cr detection in K. Solid lines are shown for ±5%, ±10%, and ±20% errors. 



5BK3er sample and 88.7% of 5B5er sample. 

We exclude photo- z measurements that are unreliable 
at the 10% level based on Az^g. This leaves us with a to- 
tal sample of 53217 galaxies. Among this accurate photo- 
z sample, 21768 galaxies are detected at K above 3a 
while the rest were not. The z p hot distribution for both 
catalogs is illustrated in Figure 6 where 27574 sources 
are identified to be at z p hot = 1-3 while 21564 are lo- 
cated below z = 1 and 4079 are above z = 3. The spike 
in the distribution at z p hot ~ 1-8 is an artifact of poorer 
sensitivity between 1 (im and 2 /zm that would capture 
the Balmcr/4000 A break, and is less significant since 
the photo- z accuracy is lower (see Figure 4). This fea- 
ture does not effect census results, since the spike is not 
located near either edge of the redshift window. 

Because of the complicated selection for the photo-z 
samples, one concern is whether derived photo-z's rep- 
resent the entire population of galaxies. The photo-z 



samples span the full range of galaxy colors observed in 
the SDF. Of course, galaxies of lower mass and luminos- 
ity are systematically missed at higher redshifts, due to 
our flux sensitivity limits. 

4.2. SED Modeling 

To determine physical properties of our galaxies, the 
Fitting and Assessment of Synthetic Templates (FAST; 
Kriek et al. 2009) code is used to model the SED. The 
spectral synthesis models are generated from the Bruzual 
& Chariot (2003) code 5 , and consist of a star forma- 
tion history that follows an exponential decay (i.e., a 
t model). For simplicity, we choose log(r/yr) = 8.0, 
9.0, and 10.0. These t values were selected to include 

5 Maraston (2005) models were also considered, but results do 
not differ significantly from those using the Bruzual & Chariot 
(2003) models. 
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Fig. 4. — Illustration of the z p hot uncertainties for the combined 
Zphot sample. The x-axis shows the best-fitting z p hot while the 
y-axis shows Azqs, the half- width of z p i lot distribution at the 68% 
confidence level normalized by l+z p hot- The horizontal lines corre- 
spond to the adopted 10% accuracy per (l + z p hot) used to exclude 
estimates. Black points denote sources in the 



'phot 



less accurate 
5BK3cr catalog. 

a bursty, intermediate, and roughly constant star for- 
mation history. In addition to the star formation his- 
tory, the grid of models spans a range of log(age/yr), 
between 7.0 and up to the age of the universe at a given 
redshift (in increments of 0.1 dex), and dust extinction 
with Ay = 0.0-3.0 mag with 0.1 mag increments. FAST 
adopts the Calzetti et al. (2000) dust extinction law, and 
we only use the solar metallicity models and assume a 
Salpeter (1955) initial mass function (IMF). 

The inputs to the SED fitting are the z p hot from EAZY, 
the photometric measurements, and the total system 
throughput for each filter (the same as those used in 
EAZY). The outputs determined from \ 2 minimization 
are the stellar masses, stellar ages, dust reddening (Ay), 
star formation rates (SFRs), and r values. Narrow- band 
filter measurements were excluded in the SED modeling 
(thus up to 17 bands), since FAST has yet to include 
emission lines. 

In the process of modeling the SEDs, we found a 
problem with the inclusion of IRAC data for 20% of 
our galaxies. First, the UV extinction-corrected 6 SFRs 
agree well with those derived from SED modeling with 
15 bands (excluding IRAC bands). However, when the 
IRAC bands were included, the same comparison of SFRs 
yielded a factor of two higher SFRs from SED modeling. 
We have not fully understood the cause(s) of this prob- 
lem, so results using information from the SED modeling 
are determined with 15 bands instead of 17 bands. 

5. PHOTOMETRIC TECHNIQUES TO IDENTIFY Z = 1-3 
GALAXIES 

The following two-color methods are used to select 
galaxies in the redshift desert: BX, BM, BzK (star- 
forming and passive), and LBG. The sizes of these galaxy 
population samples are given in Table 4. Since the J- 

6 We use the derived E(B — V) from SED modeling and assume 
a Calzetti et al. (2000) extinction law. 



band data are shallower than the K, the SDF DRG sam- 
ple alone is -~200 in size, so we do not include this in the 
census. We refer readers to Lane et al. (2007, hereafter 
L07) where a 0.5 deg 2 NIR survey was conducted and 
a large sample of bright (K < 21 AB) DRGs was ob- 
tained and compared against other NIR-selected galax- 
ies. 

5.1. BX/BM Selection 

Steidel et al. (2004, hereafter S04) defined BX and BM 
galaxies by two polygons in the U n — G/G—7Z color space. 
For BX, it is G-TZ > -0.2, U n -G > G-K+0.2, G-TZ < 
0.2([/„ - G) + 0.4, and U n - G < G - TZ + 1.0. The BM 
selection consists of G-TZ > -0.2, U n -G > G-71-0.1, 
G-TZ< 0.2(U n -G) + 0.4, and U n -G <G -TZ + 0.2. 
However, the U n GlZ system differs from existing filters 
of similar wavelengths for the SDF, so a transforma- 
tion is necessary between UBVRci' and U n GTZ. The 
transformation is derived by using Bruzual & Chariot 
(2003) spectral synthesis models of star-forming galax- 
ies. 7 Mock galaxies are produced from a grid of synthetic 
spectra between z = 0.85 and z = 3.8 and reddened with 
the Calzetti et al. (2000) law with E(B - V) = 0.0-0.4. 
We also include neutral hydrogen intergalactic medium 
(IGM) absorption following Madau (1995). Then mag- 
nitudes in the U n , G, 1Z, U, B, V, Rc, and i' filters 
are determined for these 540 mock galaxies. 8 Among the 
artificial galaxies that meet the original BX or BM selec- 
tion, least-squares fitting is used to represent the G band 
(7Z band) with a combination of B and V (Rc and i') 
such that U n - G = U - BV and G - 1Z = BV - R c i', 
where 



BV= -2.5 log 



Rci' 



-2.5 log 



xi/b + (1 - xi)fv 
3630 /iJy 

xi!r + (1 - X2)fj 

3630 ^Jy 



, and (2) 
(3) 



Here, fx is the flux density per unit frequency (erg s _1 
enr 2 Hz -1 ) in band "X" An illustration of this tech- 
nique is shown in Figure 7 for BX galaxies. 

While different values of x\ and x 2 can be adopted for 
the BX ( Xl = 0.314 and x 2 = 0.207) and BM ( Xl = 0.491 
and X2 = 0.790) selections, the net result is that some 
BX galaxies would also meet the BM selection and vice 
versa. However, the BX and BM galaxies occupy distinct 
regions in the default U n GTZ plane, so only one transfor- 
mation set can be used to identify the two galaxy pop- 
ulations simultaneously. As a compromise, we decided 
to adopt the BX x\-x 2 transformation, since it occupies 
the region between the BM and [/-dropout techniques. 
With this approach, the BX selection parallelogram re- 
mains roughly the same, while the BM and [/-dropout 
selection regions are modified. The only modification 
that we make to the BX selection is to move the lower 
bound up by 0.03 mag in U — BV: 



BV - R c i'> -0.2, 

U - BV>BV - R c i' 



0.23, 



(4) 
(5) 



7 This model is identical that used in L09 and is similar to that 
used by Steidel et al. (1999). 

8 The UnGTZ filter profiles were provided by A. E. Shaplcy. 
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Fig. 5. — Summary of the number of band detections for sources in the 5BK3<r (left) and 5B5ct (right) samples. Optical and NIR bands 
are included for a total of 15 bands. We have detections in at least 10 wave bands for 99% of 5BK3<r sample and 88.7% of 5B5cr sample. 



TABLE 4 

Summary of Photometric Selection Samples 



Method 


Section 


N 


^phot 


^phot 


2 phot< 1-0 


2phot< 0.5 


BX R c < 25.5 


5.1 


5585 


5070 


2.263±0.39 


1282 (25.3%) 


1235 (24.4%) 


BX R c < 26.0 


5.1 


8004 


6887 


2.240±0.40 


1395 (20.3%) 


1326 (19.3%) 


BM R c < 25.5 


5.1 


6408 


6097 


1.631±0.33 


943 (15.5%) 


778 (12.8%) 


BM R c < 26.0 


5.1 


8584 


7771 


1.640±0.35 


1112 (14.3%) 


862 (11.1%) 


sBzK 


5.2 


8715 


7866 


1.939±0.56 


431 ( 5.5%) 


179 ( 2.3%) 


pBzK 


5.2 


214 


208 


1.663±0.51 


1 ( 0.5%) 


( 0.0%) 


z ~ 3 LBG R c < 25.5 


5.3.1 


2375 


2060 


2.913±0.49 


307 (14.9%) 


277 (13.4%) 


2 ~ 3 LBG R c < 26.0 


5.3.1 


3823 


3265 


2.862±0.46 


348 (10.7%) 


310 ( 9.5%) 


z ~ 2 LBG V < 25.4 


5.3.2 


7878 


7663 


1.816±0.47 


619 ( 8.1%) 


238 ( 3.1%) 


z ~ 2 LBG V < 26.0 


5.3.2 


9745 


9397 


1.823±0.45 


690 ( 7.3%) 


259 ( 2.8%) 


z ~ 1 LBG 


5.3.3 


1042 


1039 


0.834±0.15 


Total (Shallow) a < b 




32217 


30003 




3583 


2707 


Unique (Shallow) a ' b ' d 




20687 


18830 




4120 


2481 


Total (Faint) a ' c 




40127 


36433 




3977 


2936 


Unique (Faint) a > c < d 




26411 


23279 




4482 


2688 



The combination of BX, BM, [/-dropout (z ~ 3), sBzK, pBzK, and TVC/V-dropout (z ~ 2) samples. 

2 BXs, BMs, and [/-dropouts selected with Rc < 25.5, and NUV- drop outs selected with V < 25.4. 

3 BXs, BMs, and [/-dropouts selected with R c < 26.0, and A/"[/K-dropouts selected with V < 26.0. 

^ "Unique" refers to the sample which accounts for overlap between the photometric selection samples. 



BV - R c i' <0-7Q(U - BV) + 0.280, and (11) 
U- BV>BV - R c i' + 0-23. (12) 

The U — BV and BV - R c i' colors for SDF galaxies 
are shown in Figure 9 for the selection of BX and BM 
galaxies. Following Stcidel ct al. (2003), we defined a 
non-detection in the U band when its flux falls below the 
lrj limit. 

In total, 6280 BX (5585 after stellar removal) and 6776 
BM (6408 after stellar removal) galaxies are identified 
with Rc < 25.5. In Section 6.2, we also discuss the BX 
and BM samples with Rc < 26.0, which consist of 8699 
BX (8004 after stellar removal) and 8952 BM galaxies 
(8584 after stellar removal). In Figure 10, the photomet- 
ric redshift distributions of the BX and BM samples are 
illustrated. Also overlaid are the spectroscopic redshift 
distributions from S04 and Reddy et al. (2008), which 
demonstrates the accuracy of EAZY at z ~ 2 and shows 
that the filter transformations did indeed approximate 
the original BX and BM definitions. Further evidence 
that the transformation was done correctly is the sur- 




BV - R c i' <0.2(U - BV) +0.4, and (6) 
U- BV<BV - R c i' + 1.0. (7) 

It should be noted that the BM transformation did 
change significantly when adopting the BX photometry 
transformation (sec Figure 8). It now spans a wider range 
in BV — Rci' color, so more z < 1 sources could poten- 
tially contaminate our BM sample at the cost of identify- 
ing z > 1 BM galaxies. As a compromise, we choose not 
to extend the selection region to encompass the lowest- 
redshift BM galaxies, in order to exclude some low-z in- 
terlopers. This may sound alarming, but it is expected, 
since the original BM selection spanned a small range 
of G — R color for a given U n — G color. Thus scatter 
from photometric uncertainties causes missed faint z > 1 
galaxies and the inadvertent inclusion of low-z interlop- 
ers. The final selection that we use for BM galaxies is 

BV - R c i' > -0.2, (8) 
U - BV>-0.1, (9) 
BV - R c i' < 0.382(C/ - BV) + 0.853, (10) 
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Fig. 6 . — Zphot distribution for the combination (solid line) of 
the 5BK3cr and 5B5<r samples and only the 5BK3cr sample (shaded 
region). Sources classified as stars were removed and objects with 
Azas> 0.1 were excluded. 

face density of BXs, which is shown in Figure 11. It is 
compared with measurements reported by Reddy ct al. 
(2008) and shows good agreement over a range of —3 
mag. 

Although not illustrated in Figure 10, a subset of the 
BX and BM samples have z pho t = 0.1-0.5. The BX 
(BM) i? c < 25.5 sample has 1235 (778) with z phot < 0.5 
corresponding to 24.4% (12.8%). These numbers de- 
crease to 19.3% (11.1%) when considering the sample 
with Rq < 26.0. These galaxies have accurate z p hot's 
(Az68 < 4%), so it cannot be argued that these galaxies 
are at z > 1. Furthermore, we inspected the B — z' and 
z 1 — K colors and find that these objects occupy a region 
below the star-forming BzK selection (see Section 5.2), 
which provides further evidence that they are at z < 1.3 
(i.e., the Balmer/4000 A break is observed in the opti- 
cal) . The full z spoc distributions for BX and BM galaxies 
were not provided in S04, but the z spcc < 1.0 interloper 
statistics were stated to be 17.1% and 5.8% for BX and 
BM (1Z < 25.5) galaxies, respectively. G07 found similar 
results to ours for their BX sample. 

5.2. BzK Selection 

Daddi et al. (2004) defined star-forming BzK (sBzK) 
galaxies as those with BzK = (z — K) — (B — z) > —0.2 
with detection required in B and K. This corresponds to 
the solid line in Figure 12, which has 8715 star-forming 
BzK (hereafter sBzK) galaxies falling above it. In addi- 
tion, the passive BzK (pBzK) selection, which requires 
BzK < -0.2 and z - K > 2.5, yields 228 galaxies (214 
after stellar removal) with z = 1-1.5. The photo-z's of 
the pBzK and sBzK populations arc shown in Figure 13. 
They show good agreement with photo-z's derived for 
BzK galaxies in the COSMOS field (McCracken et al. 
2010). The COSMOS BzK sample probed K < 23.0 
while our sample reaches K ~ 24.0. We find better 
agreement with COSMOS at the z p hot = 1-7-3.0 range 
when limiting our sample to K < 23.0. The COSMOS 



and SDF samples of BzK galaxies independently confirm 
that sBzK galaxies are at z w 1.5-2.5 and pBzK galaxies 
are at z w 1.0-1.5. 

In Figure 14, we plot the surface density of the SDF 
sBzK galaxies and previous sBzK surveys (Reddy et al. 
2005; Kong et al. 2006; Lane et al. 2007; McCracken et 
al. 2010). We find that our measured surface density 
for sBzK is consistent with those reported in McCracken 
ct al. (2010) down to their if-band depth of -23.0 AB. 
Kong et al. (2006) show a higher surface density but is 
consistent with the photometric scatter. L07 also find a 
higher surface density, although McCracken ct al. (2010) 
pointed out that this is likely due to a photometric cali- 
bration issue. Finally, the surface density of R05's sBzKs 
is consistent with ours for K < 23 AB, and then is lower 
than ours by a factor of at most two. We will discuss 
this discrepancy at faint K in Section 6.3. 

5.3. LBG Selection 

The above techniques were designed to select z = 1- 
3 galaxies when the Lyman break technique is not de- 
tectable with optical photometry. This has now changed 
since L09 used ultraviolet imaging for the first survey 
of z ~ 2 LBGs in the SDF. To fully span the z = 1 
3 era with the Lyman break technique, we select U-, 
NUV-, and ft/F-dropouts. We will interchangeably use 
U/NUV/FUV-dmpouts and z - 3/2/1 LBGs to refer to 
the Lyman-limit break selected sources. 

5.3.1. U -dropouts 

Using the UBVRci' to U n GlZ transformation de- 
scribed in Section 5.1, we identify [/-dropouts in a sim- 
ilar manner as Steidel et al. (2003): U - BV > -0.1, 
BV - R c i' < 0.95, and U - BV > BV - R c i' + 1.0. 
A modification is made to the right edge by shifting it 
0.15 mag bluer in BV — Rci' ■ This selection spans the 
region above the BX selection window, as illustrated in 
Figure 15. We find 2578 [/-dropouts with R c < 25.5 
(2375 after stellar removal). The surface density for our 
U -dropouts (see Figure 11) is similar to those reported by 
Reddy et al. (2008). Finally, we compare the z p hot 

for U- 

dropouts against the Reddy et al. (2008) sample in Fig- 
ure 10 and find good agreement, which further supports 
the idea that our selection is probing the same galaxies 
identified by Steidel et al. 

5.3.2. NUV-dropouts 

The selection of ^[/U-dropouts is discussed in L09. 
The definition is: NUV - B > 1.75, B — V < 0.50, 
and NUV -B> 2.4(5 - V) + 1.15. In defining the se- 
lection criteria, we used a spectral synthesis model that 
is similar to that used for selecting [/-dropouts. The 
range of dust reddening probed by the above color cri- 
teria, E(B — V) = 0.0-0.4, is fairly similar to those se- 
lected by Steidel et al. should galaxies that span this 
full range exists. Thus, the main difference between the 
two populations of LBGs is the redshift sensitivity. This 
statement, of course, assumes that there is no evolution 
in the properties (e.g., dust) of LBGs between z ~ 3 and 
z - 2. 

There are two corrections that we make to the previous 
selection. First, a source is considered undetected in the 
NUV when it is below the la threshold. Previously, we 
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Fig. 7. — BX transformation for U n GlZ to U BV Rci' ■ (Left) The originally defined selection of BX galaxies in the U n GlZ space. (Right) 
The U — BV and BV — Rci' colors for these BX galaxies, as determined by Equations (2) and (3) with x\ = 0.314 and 1E2 = 0.207. The 
solid black lines refer to the final BX selection (see Equations (4)-(7)) while the red dashed lines illustrate the S04 selection. These plots 
indicate a proper transformation has been made to select all BX galaxies. These transformation is further discussed in Section 5.1. (A 
color version of this figure is available in the online journal.) 
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Fig. 8. — Similar to Figure 7, but illustrated here for BM galaxies. These plots indicate a reasonable transformation for BM galaxies when 
using the BX transformation except for z ~ 1 BM galaxies. We adopt a different selection for BM galaxies, as shown by the solid black 
lines (see Equations (8)-(12)) compared to the S04 selection illustrated by the red dashed lines. Further discussion of the transformation 
is available in Section 5.1. (A color version of this figure is available in the online journal.) 



adopted 3<r to be conservative; however, we decided to 
be consistent with Steidcl ct al. (2003) for the final selec- 
tion of all Lyman-limit galaxies. This will shift points in 
Figure 9 of L09 by 1.2 mag redder in NUV — B. Second, 
we failed to apply an aperture correction for undetected 
sources in L09. The aperture correction that we adopted 
previously was 1.81, and our latest calculation, which we 
now use, indicates 1.85, which leaves a 0.6 mag offset in 
the NUV — B colors for faint galaxies. 

By adopting a la threshold in the NUV, the z ~ 2 UV 
luminosity function discussed in L09 changes. To first 
order, the completeness corrections were underestimated 
at z = 2-2.5, since some of the simulated galaxies that 
fell below our NUV — B = 1.75 mag will now enter the 
z ~ 2 LBG selection region. However, by adopting a 
lower threshold, more sources will also enter our selection 
region that were bluer than NUV — B = 1.75 mag. We 
note that the comparison between BX/BM and z ~ 2 
LBG is more robust in this paper than in L09 since a 
direct comparison is made using the same data. 

The NUV - B and B — V colors for V < 25 A sources 



are shown in Figure 16, and 8574 ./Vf/F-dropouts (7878 
after stellar removal) are identified. Recall that the NUV 
fluxes were obtained by PSF-fitting bright sources in or- 
der to obtain reliable fluxes for faint, confused sources. 
If this step was not performed, the A/f/V-dropout sam- 
ple would be ~30% smaller. This is similar to the frac- 
tion of missed mock LBG galaxies in the Monte Carlo 
simulation of L09. In Section 6.2, we also discuss the 
AffJK-dropout sample with V < 26.0, which consists of 
10441 sources (9745 after stellar removal). The z p hot dis- 
tribution is shown in Figure 16, and indicates that the 
selection mostly identifies galaxies at z = 1.5-2.5, with 
25.4% contamination from z < 1.5 interlopers. 

5.3.3. FUV-dropouts 

Studies (e.g., Burgarella et al. 2007) have used GALEX 
FUV imaging to select z ~ 1 LBGs. The selection of 
.Ft/V-dropouts is currently best described in conference 
proceedings. The selection which has been adopted is 
FUV - NUV > 2.0, NUV - U < 1.4, and FUV ~ 
NUV > 1.05(NUV - U) + 2.52. The FUV - NUV 
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Fig. 9.— U — BV and BV - R c i' colors to select BX and BM 
galaxies. Triangle points are sources that are undetected in U at 
lcr. Black lines indicate the selection boxes for BX (top) and BM 
(bottom). Selection is performed for Rc < 25.5 mag. 
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Fig. 10. — Distribution of z phot 's for the UV-selected BX, BM, 
and [/-dropout (solid lines) compared to the spectroscopic samples 
of S04 or Reddy ct al. (2008). Galaxies were selected down to 
Rq = 25.5 mag. 

and NUV — U colors are shown in Figure 17. We iden- 
tified 1067 FC/U-dropouts (1042 after stellar removal) 
with U < 25.0 mag. We find that the z p hot distribution 
is skewed towards z ~ 0.7-0.8 (see Figure 17). While 
Ff7V-dropouts have been claimed to be at z ~ 1, the Ly- 
man continuum breaks occurs at the FUV filter center 
at z = 0.68, reducing the flux by half. This peak is not 
surprising given (1) the current UV sensitivity and (2) 
the majority of z ~ 1 galaxies have U > 25 mag. 

5.4. Summary 



We used color selection to identify 5585 BXs (8004 
R c < 26.0), 6408 BMs (8584 R c < 26.0), 8715 sBzKs, 
214 pBzKs, 2375 [/-dropouts (3823 Rc < 26.0), 7878 
AW-dropouts (9745 V < 26.0), and 1042 F£/U-dropouts 
for a total of 32217 (adopting a shallower limit on Rc 
or V) or 40127 (adopting a fainter limit on Rc or V) 
galaxies potentially at z « 1-3. However, since there is 
overlap between the different color-selected samples, the 
non-redundant sample contains 20687 (18830 with z p hot) 
and 26411 (23279 with z p hot) galaxies in the Shallow and 
Faint samples, respectively. 

6. RESULTS 

In this section we begin by discussing the properties 
of galaxies selected using a given photometric technique. 
Then we compare the different color selections to dis- 
cuss sample overlap and the selection bias associated 
with each technique. We then compare the color-selected 
samples against photo-z's to estimate the completeness 
of using these techniques in acquiring a census of star- 
forming galaxies at z = 1-3. 9 

6.1. The Physical Properties of z = 1-5 Galaxies 

Using the color selection samples with z p hot and infor- 
mation from SED modeling, we illustrate in Figure 18 
the efficiency of these photometric techniques in terms 
of the relative fraction of the total number. We find that 
the BX and BM methods identify 60%-80% of galax- 
ies at the peak of their z p hot distributions (i.e., where 
they are most sensitive). The [/-dropout technique iden- 
tifies 80% of galaxies at z ~ 3. However, this is op- 
timistic since our "Total" sample is missing low-mass 
z ~ 3 red/dusty galaxies should they exist. Likewise, 
it appears that the iV[/U-dropout technique has similar 
efficiency as the BX/BM method with a redshift coverage 
that spans the redshift of BX and BM galaxies. LBGs at 
z ~ 1 would require UV imaging at «1800 A (between 
the FUV and NUV band). Finally, the star- forming BzK 
selection probes a wider range in redshift: capturing at 
least 40% of galaxies at z > 3 and as much as ~70% at 
z ~ 1.5. 

Figure 18 also illustrates that the very different meth- 
ods to identify galaxies at this epoch are complemen- 
tary, and that no single method identifies all galaxies. 
The union of the different color selection techniques does 
yield a comprehensive galaxy census of z = 1-3. 

To further understand the selection bias of these tech- 
niques, we illustrate in Figure 19 the distribution of dust 
reddening, stellar mass, SFR, and stellar ages for the BX, 
BM, sBzK, pBzK, and NUV-diopout samples. These 
distributions are normalized to the census of these tech- 
niques. The majority of galaxies with stellar masses 
above 6 x 10 9 M Q are acquired with the sBzK tech- 
nique. UV selection methods are more efficient at se- 
lecting lower-mass galaxies. The BzK technique is also 
able to identify the oldest galaxies and those with high 
dust extinction. UV techniques find young galaxies with 
little dust extinction and typical SFRs of 1-10 Mq yr -1 . 
The z ~ 2 LBG population appears to span a wider range 
of dust reddening compared to the BX/BM method. We 
attribute this to the (1) larger range in rest-UV colors 

9 We exclude the FC/V-dropouts since the sample overlap with 
others is small and the z p hot distribution peaks at z < 1. 
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Fig. 11. — Observed surface density of BX galaxies (left) and [/-dropouts (right) shown as filled (exclude stars) and open (include stars) 
squares. There is good agreement with the filled circles from Reddy et al. (2008). 
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Fig. 12. — B — z' and z' — K colors for ~25,000 sources detected 
in the K band above 3cr. Star-forming BzKs are identified as those 
above the slanted line (BzK = —0.2) while the selection of passive 
BzKs is through a combination of the horizontal and the slanted 
lines. Triangle points indicate A"-band sources that are undetected 
in the B band at 3<r (lower limit on the B — z' color). 

(B — V < 0.5 mag spans a wider G — 1Z color compared 
to the G — TZ color of BX/BM galaxies) and (2) an abil- 
ity to select massive galaxies. We discuss this further in 
Section 6.2.1. Figure 19 also shows that the most ac- 
tively star-forming galaxies 10 will be identified using ei- 
ther the sBzK or A^t/V-dropout method. However, 66% 
of the color-selected z p hot = 1-3 census consists of galax- 
ies with SFRs between a few and 30 M Q yr _1 , and 70% 
of this fraction is captured using the BX/BM method 
with Rq < 25.5 mag. 

6.2. A Direct Comparison of Color Selection Techniques 



10 The high SFRs are due to large dust extinction corrections. 



1 UDS is short for the Ultra Deep Survey, which is part of the UKIRT 
Infrared Deep Sky Survey. 

^ The SDF z Bpcc is mostly limited to low-2, so nothing is reported 
here. 

^ The number here is merged with the above number for BX. 
^ EROs (extremely red objects) are often distinguished in the litera- 
ture from DRGs. This selection consists of using an i? — K > 5.3 mag 



Previous studies (R05; Q07; G07; L07) have also in- 
vestigated the completeness of photometric selection and 
the selection bias in the GOODS-N, GOODS-S/CDF-S, 
MUSYC, and UKIRT fields. However, as illustrated in 
Table 5, each study has its own advantages and disadvan- 
tages in terms of area coverage, color selection used, and 
redshift accuracy. Since these surveys could not iden- 
tify z ~ 2 LBGs, comparisons were made between the 
BX/BM/cT-dropout methods and the BzK /DRG meth- 
ods given, and it was assumed that BX/BM's arc lowcr- 
rcdshift analogs of z ~ 3 LBGs. We discuss the overlap 
between BX/BM and z - 2 LBGs in Section 6.2.1. In 
previous studies, the term "LBG" is used to refer to the 
selection of [/-dropouts. We make a clear distinction be- 
tween the two. 

A summary of the color selection sample overlap is pro- 
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Fig. 13 . — Zphot distribution for 7866 of 8715 sBzK galaxies (left) and 208 of 228 pBzK galaxies (right) indicated by the solid lines. 
The dotted line represents sBzK galaxies with K < 23 mag to compare with the photo- z distribution of sBzK galaxies from the COSMOS 
survey (McCracken et al. 2010), which are overlaid as dashed lines. The COSMOS distributions have been normalized to the peak in the 
SDF K < 23 mag z p hot distributions. 
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Fig. 14. — Surface density of sBzK galaxies in the SDF (filled cir- 
cles) compared to other surveys: McCracken et al. (2010, squares), 
L07 (triangles), Kong et al. (2006, asterisks), and R05 (diamonds). 
Red error bars illustrate the minimum and maximum surface den- 
sity for eight independent 8f5x8.'5 fields in the SDF. (A color ver- 
sion of this figure is available in the online journal.) 
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Fig. 15. — U — BV and BV — Rci' colors to select [/-dropout 
galaxies. Triangle points are sources that are undetected in U at 
lcr. The selection boxes for BXs, BMs, and {/-dropouts are shown. 

with such restrictions in Section 6.2.1, to compare with 
previous estimates. 



vided in Table 6, along with the average z p hot, E(B — V), 
stellar age, stellar mass, and SFR. We have created a 6- 
digit binary index to refer to the sample overlap. Each 
digit refers to a specific galaxy population in the fol- 
lowing order: pBzK, sBzK, BM, BX, [/-dropout, and 
./VE/V-dropout. So for example, galaxies that simultane- 
ously satisfy the color criteria of sBzK, BX, and NUV- 
dropout are denoted as "010101". The diagram in Fig- 
ure 20 illustrates the sample overlap. The BX, BM, and 
[/-dropouts do not overlap, given the mutually exclusive 
color selection. This is also the case for the passive and 
star-forming BzK-selected galaxies. Note that we have 
not placed any restrictions on the optical (if-band) mag- 
nitudes for NIR-selected (UV-selected) galaxies for this 
table and this overlap diagram. We discuss the overlap 



6.2.1. Sample Overlap: Color-Color Plots 

To further understand and illustrate the selection ef- 
fects of these techniques, we plot where each galaxy pop- 
ulation/technique (e.g., BX) falls in the color-color se- 
lection of another (e.g., BzK). We first discuss the BzK 
selection, followed by the U n GlZ methods, and finally 
the selection that we have developed to identify z ~ 2 
LBGs. Note that the numbers below have excluded in- 
terlopers with Zphot < 0.5 (unless otherwise indicated). 
When computing the overlap fraction and providing the 
two-color distribution of sources, we place limits on ei- 
ther the V, Rq, or K, as we did when we generated the 
high-z photometric samples. 

A detailed summary of the statistics reported in this 
section is provided in Tables 7-9. 
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Fig. 16. — (Left) The selection of iVfTV-dropouts with NUV — B and B — V colors with V < 25.4. Triangle points indicate sources with 
lcr limits in the NUV band. (Right) The z p hot distribution for most ./VCTV-dropouts. 
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Near-IR selection. The B — z' and z' — K colors for 
BXs, BMs, ./V[/U-dropouts, and {/-dropouts arc shown 
in Figure 21. When the BX, BM, [/-dropout, and NUV- 
dropout samples (1) are restricted to the region of sensi- 
tive K photometry, (2) exclude sources with z pn ot < 0.5, 
and (3) only consider sources detected above 3er in K, 
we find that 96.6% of BXs, 76.6% of BMs, 67.8% of U- 
dropouts, and 75.3% of iVW-dropouts meet the sBzK 
selection, which indicates that sBzK method is capable 
of selecting almost all high-z UV-selected galaxies hav- 
ing K < 24 mag. Table 7 provides more information 
regarding these statistics. We note that if we did not 
exclude objects with z p hot < 0.5, the overlap fraction 
would be 75.9% (BX), 68.5% (BM), 59.6% ([/-dropout), 
and 73.5% (AW-dropout). More than two-thirds of U- 
dropouts with if-band detections are also identified as 
sBzK. The redshift where the Balmer/4000 A break en- 
ters the if-band window is z w 4, so extremely deep K- 
band imaging (K = 25-26) can yield mass-limited sam- 
ples of z < 4 galaxies with stellar masses above ~ 10 10 



Mq. These high-z sBzKs were also seen by G07, but 
in smaller numbers. Unfortunately, they did not report 
a sBzK-[/-dropout overlap fraction to compare to ours. 
G07 did find that > 86% of z > 1.4 BX/BM galaxies are 
sBzK galaxies, which is consistent with the above num- 
ber. Q07 determined that 80% of BX/BM/ [/-dropout 
galaxies (regardless of z p hot) arc sBzK-selcctcd, which 
agrees with our number without any z p hot constraints. 
Likewise, similar numbers were reported by R05. 

The majority of sources that are missed by the sBzK 
method are BM and ./V[/F-dropout galaxies. We find that 
these non-sBzK galaxies have characteristic Zphot ~ 1.9, 
E(B - V) w 0.1, stellar masses of - 10 10 M , SFRs of 
~10 Mq yr _1 , and stellar ages of 3 x 10 8 yr. It is not 
surprising that the BzK technique misses blue galaxies, 
since it is more sensitive to redder galaxies. 

None of the BX, BM, [/-dropout, and iV[/U-dropout 
sources are classified as a passive BzK, indicating that 
pBzK galaxies represent a completely separate popula- 
tion. This was also seen by Q07. We further discuss this 
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14±0.10 


10 


07±0.42 


1.33±0.44 


8.53±0.46 


010011 


77 


2 


50±0.35 





13±0.11 


10 


88±1.28 


1.65±0.75 


8.74±0.48 


000001 


1266 


1 


23±0.40 





27±0.16 


9 


78±0.82 


1.21±0.82 


8.19±0.70 


001001 


1470 


1 


73±0.26 





15±0.09 


9 


58±0.80 


1.14±0.54 


8.23±0.80 


000101 


1043 


2 


29±0.29 





10±0.08 


9 


50±0.81 


1.15±0.45 


8.15±0.78 


000011 


95 


2 


59±0.25 





08±0.11 


9 


82±1.46 


1.38±1.12 


8.14±0.63 


001000 


1622 


1 


47±0.40 





12±0.09 


9 


17±0.62 


0.71±0.48 


8.26±0.66 


000100 


908 


2 


30±0.53 





10±0.10 


9 


43±0.77 


0.86±0.50 


8.22±0.70 


000010 


1167 


2 


99±0.51 





07±0.09 


9 


75±0.69 


1.05±0.59 


8.33±0.54 


Note. 


— Average properties arc reported. 


The first 


column is a 6-digit binary 



flag with each digit referring to a specific galaxy population in the order of pBzK, 
sBzK, BM, BX, (/-dropout, and TVL^V-dropout. See Section 6.2 for further description. 
Sources with £ p hot < 0.5 arc excluded. 



N = 200, each bin 




1.0 1.5 2.0 2.5 3.0 
Photometric Redshift 



3.5 



Fig. 18. — Fraction of the Shallow sample versus redshift for the 
different color selection techniques: AfC/V-dropout (solid; black) 
sBzK (dashed; red), pBzK (dotted; orange), BX (dot-dashed 
green), BM (dot-dot-dashed; blue), and [/-dropout (long-dashed 
purple). We illustrate (in a given redshift bin) the fraction in terms 
of the total number of galaxies. (A color version of this figure is 
available in the online journal.) 
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Fig. 19. — Physical properties of z = 1-3 galaxies. The 
distribution of A v (top left), log(SFR[Af Q yr - *]) (top right), 
log(Af s tar[AfQ]) (bottom left), and log(age[yr]) (bottom right) for 
the different population. Color and line-style conventions follow 
those in Figure 18. The [/-dropouts are not shown since most of 
them are at z p hot ^ 3.0. (A color version of this figure is available 
in the online journal.) 



TABLE 7 

Summary of Sample Overlap: NIR Selection 



Type 


FoV Limited 


K > 3a 


Total 


BX 


3493 


2065 


1995 (96.6%) 


BM 


4670 


2955 


2265 (76.6%) 


z~3 LBG 


1701 


824 


559 (67.8%) 


z~2 LBG 


6249 


4784 


3602 (75.3%) 



Note. — Sources with z p hot< 0.5 have already been 
excluded. 



TABLE 8 

Summary of Sample Overlap: U n G7l Method 



Type 


R < 25. 5 a 


BX 


BM 


z~3 LBG 


Total 


sBzK 


6201 


1995 


2265 


559 


4819 (77.7%) 


z~2 LBG 


7603 


2465 


3151 


180 


5796 (76.2%) 



Sources with z p hot < 0.5 are excluded. 

in the Appendix. 

The fraction of sBzK galaxies that arc BXs, BMs, or 
[/-dropouts varies with if-band magnitude, with larger 
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BM BX z = 3 LBG 



pBzK 



100000 0213 

0212 V >= 25.4 
0162 R c >= 25.5 
of Total: 1.23% 
of pBzK: 99.53% 

z = 2 LBG 




001000 1838 

1245 K < 3ct 
0881 V >= 25.4 
of Total: 10.61% 
of BM: 32.65% 

K flag: 404/1838 




000100 1257 

1026 K < 3(7 
0817 V >= 25.4 
of Total: 7.26% 
of BX: 28.90% 

K flag: 390/1257 




000010 1440 

0995 K < 3cr 
0705 V >= 25.4 
of Total: 8.31% 
of z = 3: 68.64% 

K flag: 366/1440 




000001 1313 

0216 K < 3ct 
0036 R c >= 25.5 
of Total: 7.58% 

of z = 2: 17.19% 
K flag: 337/1313 


001001 1527 

0809 K < 3(7 

of Total: 8.82% 
of BM: 27.12% 
of z = 2: 19.99% 
K flag: 556/1527 




000101 1098 

0724 K < 3(7 

of Total: 6.34% 
of BX: 25.24% 
of z = 2: 14.37% 
K flag: 467/1098 




00001 1 0099 

0068 K < 3(7 

of Total: 0.57% 
of z = 3: 4.72% 
of z = 2: 1 .30% 
K flag: 31 /99 




010001 0530 

of Total: 3.06% 

of z = 2: 6.94% 
of sBzK: 6.21% 


011001 1624 

of Total: 9.38% 
of BM: 28.85% 
of z = 2: 21 .26% 
of sBzK: "9.03% 




010101 1367 

of Total: 7.89% 
of BX: 31 .43% 
of z = 2: 1 7.89% 
of sBzK: 1 6.0" % 




01001 1 0081 

of Total: 0.47% 
of z = 3: 3.86% 
of z = 2: 1 .06% 
of sBzK: 0.95% 


010000 3187 

2843 V >= 25.4 
2334 R c >= 25.5 
of Total: 18.40% 

of sBzK: 37.34% 


01 1000 0641 

0218 V >= 25.4 

of Total: 3.70% 
of BM: 1 1 .39% 
of sBzK: 7.51% 




010100 0628 

0239 V >= 25.4 

of Total: 3.63% 
of BX: 14.44% 
of sBzK: 7.36% 




010010 0478 

0171 V >= 25.4 

of Total: 2.76% 
of z = 3: 22.78% 
of sBzK: 5.60% 















Fig. 20. — Diagram illustrating photometric galaxy sample overlap for the Shallow sample. Redshift increases to the right. Statistics are 
provided in each box including the size of a given sample, and the fraction relative to the full and each individual sample. "K flag" refers 
to sources that fall in lower sensitivity regions, so they cannot be classified in the BzK diagram. We exclude sources with z p ^ ot < 0.5. 



TABLE 9 

Summary of Sample Overlap: A?"(_7V-dropout 
Method 



Type 


%>hot > 1.5 


V < 25.4 


Total 


BX 


3697 


2934 


2337 (79.7%) 


BM 


3499 


2886 


2512 (87.0%) 


sBzK 


6300 


4095 


3160 (77.2%) 



overlap at fainter K (see Figure 22). This is to be ex- 
pected, as the UV selection techniques probe less red- 
dened and young galaxies that are typically of lower 
masses. Our measurements agree with those of R05 at 
K « 22-23 mag and extend toward K « 24 to show a 
slightly larger fraction of overlap. However, our sBzK 
BX/BM/[/-dropout overlap fraction is noticeably lower 
for K < 22.0 mag. We find that the overlap fraction for 
a magnitude bin of 20 < K < 22 is typically 23% with lcr 
variation of 10% when dividing our SDF survey into cells 
of 8'5 x 8.'5, which are similar to the surveyed area of 
R05. Thus, the discrepancy could be explained by clus- 
tering of such galaxies and their smaller area coverage. 

Non-ionizing UV continuum selection: the U n GlZ 
method. For the selection of BX, BM, and [/-dropout 
galaxies, the BzK galaxies and z ~ 2 LBGs are plotted 
on the U — BV and BV — Rci' color space in Figure 23. 
We find that among the sBzK galaxies with Rc < 25.5 
and Zp hot > 0.5, 32.2% are BXs, 36.5% arc BMs, and 
9.0% are [/-dropouts, for a total overlap of 77.7%. De- 
tailed information of the overlap is provided in Table 8. 



Q07 determined ^65% for this fraction. To understand 
this small yet noticeable difference, we point out that our 
sample probes much fainter K magnitudes than they did, 
by «1 mag. Since the majority of BX/BM galaxies have 
low stellar masses, many of them arc faint in K . As a re- 
sult, our survey is more complete in mass, and is weighted 
more toward the lower luminosity galaxies. We see this 
from the slightly higher overlap fraction in Figure 22 for 
A'ab > 23.0 mag. Likewise, among iVf/V-dropouts with 
R c < 25.5 and z phot > 0.5, 32.4% are BXs, 41.4% are 
BMs, and 2.4% are [/-dropouts, for a total overlap of 
76.2%. The small overlap between iV[/l^-dropouts and 
[/-dropouts is to be expected, since we intended to ex- 
clude z > 2.7 galaxies with the B — V < 0.5 mag criterion 
(see L09). 

Recall that S04 defined the BX/BM selection to 
identify z ~ 2 galaxies that are analogous to z ~ 3 
LBGs. The presence of this large overlap shows that 
this statement is mostly true. We find, however, that 
the LBG selection identifies a significant population 
of redder and more massive galaxies: at least half of 
the K < 21.5 sBzK galaxies are selected using the 
/Vf/l^-dropout method. This was previously illustrated 
in Figure 21 with the iVfJF-dropouts spanning redder 
B — z' and z' — K colors compared to the BX galaxies, 
and again in Figure 22, where the iVi/V-dropout method 
identified more K bright galaxies compared to BX/BM. 
It also explains the cluster of points at BV — Rci' = 0.5- 
1.3 mag, which is also seen for sBzK galaxies. 
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BX Galaxies 



BM Galaxies 





2 - 



tiitisP? 



o - 






-10 12 3 

B — z (AB) 

Fig. 21. — Sample overlap: B — z' and z' — K colors for BX galaxies (upper left), BM galaxies (upper right), WC/V-dropouts (lower 
left), and t/-dropout galaxies (lower right). Sources shown as black are those with > 3cr detection in K, while those in gray are sources 
undetected in K (upper limits on their z' — K color). Note that the axis scales differ from those in Figure 12. 

Ionizing UV continuum selection: the NUV-dropout 
method. We show the NUV-B and B-V colors for BX, 
BM, and BzK galaxies in Figure 24. We limit the BX, 
BM, and sBzK samples to sources brighter than V = 25.4 
and require that z p \ wt > 1.5, since the TVf/Vwill only be 
sensitive to a Lyman break above this redshift. We find 
the fractions of BX, BM, and sBzK samples that meet 
the TVf/V-dropout criteria are 79.7%, 87.0%, and 77.2%, 
respectively. More information regarding these statistics 
are given in Table 9. The high overlap strongly suggests 
that using a Lyman-limit break to select high-z galaxies 
of low and high stellar masses is a successful alternative 
to those methods. 

Summary. We find that the overlap fraction between 
UV- and NIR-selected samples (with typical restrictions 
on optical and NIR depths) of z = 1-3 galaxies is high: 
74%-98%. The sample overlap was greater when exam- 
ining the UV-selected samples in the BzK color space 
rather than vice versa. This illustrates that UV selec- 
tion techniques, which are good at identifying blue high-z 
galaxies, miss the redder galaxies, particularly the more 
massive ones. This has been seen in previous studies 
(e.g., Q07). Several comparisons between the BX/BM 



galaxies and the z ~ 2 LBGs confirm that most (80%- 
87%) BX/BM galaxies do show a strong Lyman-limit 
break. This implies that these UV-selected galaxies are 
nearly analogous to Lyman break selected galaxies. How- 
ever, the VC/V-dropout population docs span a wider 
range in dust extinction, and shows greater overlap with 
bright sBzK galaxies. 

6.3. The Star Formation Rate Density from the Census 

With a census of star-forming galaxies from the com- 
bined color selection techniques, we can determine the 
SFR density from the results of our SED modeling with 
FAST. R05 first estimated this with a census of DRGs, 



spec 



< 



BX and BM galaxies, and sBzK galaxies at 1.4 < z, 
2.6. We consider sBzK, BX and BM, and iVf/V-dropout 
galaxies, and limit our sample to 1.5 < z p hot < 2.6, since 
this is the overlap range where all three techniques work. 
While our census does not include DRGs, R05 reported 
that these galaxies contribute < 20% towards the SFR 
density via their census. 

We determine the integrated SFR density as a function 
of K- and i?c-band magnitudes for the three populations 
and the full census accounting for sample overlap. To 
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Fig. 22. — Fraction of z p hot > 1-5 sBzK galaxies that meets 
the BX (diamonds), BM (triangles), (7-dropout (squares) or NUV- 
dropout (crosses) selection. The combination of the first three 
UV selection techniques is shown as filled circles. The sample 
overlap increases towards lower X-band luminosities, which is ex- 
pected given the selection bias that UV techniques target lower 
mass galaxies. Overlaid in asterisks are estimates from R05. This 
figure also illustrates that the Lyman break technique is able to 
identify many massive z ~ 2 galaxies, which the BX/BM method 
misses. 

normalize our SFR measurements by the effective sur- 
vey volume, we assume that the completeness C(z) of 
the census at any given redshift is proportional to the 
Zphot distribution normalized to unity at the z p hot peak. 
We note that this is a lower limit on incompleteness, so 
SFR density measurements reported here can be under- 
estimated. The following equation, 



dzC(z) 



z=1.5 



dV 
dzdfl ' 



(13) 



yields volumes of 1.18 x 10 6 Mpc 3 and 1.24 x 10 6 Mpc 3 
when we consider sources with K < 24 and Rq < 25.5 
mag, respectively. We determined a total SFR density 
of 0.18±0.03 (0.18±0.03) M yr -1 Mpc- 3 for sources 
brighter than K — 24.0 mag (i?c = 25.5 mag). The 
SFR density is illustrated in Figure 25, where the range 
of allowed values is la from including uncertainties in 
the SFRs due to the z p hot accuracy and cosmic variance 
(estimated from Somerville et al. 2004). 

It appears that the SFR density down to a stellar 
mass limit is mostly contributed by sBzK galaxies. They 
are not only numerous in the sky, but also have signif- 
icant star formation that is obscured by dust. For ex- 
ample, while one-fourth of the census is composed of 
star-forming galaxies with E(B — V) > 0.25, they ac- 
count for 65% of the integrated SFR density. If a cut of 
E(B -V) = 0.4 mag (recall that the BX/BM and NUV- 
dropout techniques are sensitive up to E(B — V) ~ 0.4 
mag) is adopted, then «64% of the total SFR density is 
from "less" dusty galaxies. 

Compared to R05, we find a higher integrated SFR 
density by 67%. However, R05 was unable to select 
iVW-dropouts. Therefore, if we removed galaxies that 



uniquely meet the iVf/V^-dropout selection, then the SFR 
density for K < 24.0 mag is reduced by 10% to 0.16±0.03 
Mq yr -1 Mpc -3 , which is still ^50% higher than the re- 
ported value of R05. 

One explanation for the comparatively higher SFR 
density is the effect of ficld-to-ficld fluctuations on the 
sBzK galaxy population. In Figure 14, we compared dif- 
ferent sBzK surveys and found that the SDF sBzK counts 
agree well with those from COSMOS (McCracken et al. 
2010) down to the COSMOS survey limit of K ~ 23. 
However, R05 begin to see a relative decline at K ~ 23 
and is approximately a factor of two lower than our mea- 
surements at K oj 24. When limited to z = 1.4-2.6, 
the surface density reported by R05 is 4.93 arcmin -2 
for K < 24, which is 27% lower when compared to the 
SDF down to the same limits and redshift range (6.78 
arcmin -2 ). 

To understand the general discrepancy in surface den- 
sity, we reasoned that field-to-field fluctuations is a likely 
culprit since R05 covers 72.3 arcmin -2 . As an exercise, 
we subdivided the SDF into eight independent fields each 
with the same surveyed area of R05. We determined the 
surface density as a function of K for each field, and 
found that the variations observed can fully explain the 
discrepancy for K < 23 and about half of the differ- 
ences for K ~ 23-24. If this explanation is correct, then 
such variations can easily affect both the cumulative SFR 
density from the census of color selection techniques and 
the contribution to the census solely from sBzK galax- 
ies. For example, among the BX/BM and sBzK census, 
we find that 66% and 89% are from BX/BM and sBzK 
galaxies, while R05 report 68% and 67%. In addition, 
the smaller area coverage of R05 yields a census sample 
of ^500 galaxies in size, which is statistically less robust 
when compared to the SDF sample of ^5200 for z = 1.5- 
2.6. 

6.4. The Completeness of Color Selection Techniques 

In Figure 26, we show the redshift distribution for 
the different color selection techniques and the photo- 
z sample. For this comparison, wc consider photo-z 
sources that meet the selection depth used for color se- 
lections (i.e., 3ct detection in K 7 Rq < 25.5 mag, or 
V < 25.4 mag). This yields a total of 17566 galaxies 
with Zphot = 1.0-3.0. It is apparent from Figure 26 that 
the color selection techniques each identify some subset 
of galaxies with z p hot = 1 _ 3. We find that the combi- 
nation of different color selection techniques is able to 
identify «90% of galaxies (down to the above magnitude 
limits) with z p hot = 1.5-2.5. Since the photo-z sample 
is limited to the depths of the color selection techniques, 
the missed rate of 10% is likely a result of photometric 
scatter that drives these galaxies out of the color selec- 
tion regions. 

As expected, a noticeable («50%) fraction of galaxies 
with z p hot = 1.0-1.4 are missed. This is because (1) the 
BzK selection was designed to search at z > 1.5, (2) two- 
thirds of the BM galaxy population lie at z p hot > 1-5, 
(3) almost all BX galaxies are above z p hot > 1-5, and (4) 
the GALEX/NW filter only detects a Lyman continuum 
break at z p hot > 1.5. However, as we discussed previ- 
ously, recently available Hubble/WYCi filters are begin- 
ning to be used to select galaxies via the Lyman break 
at z ~ 1.3, and such samples may be able to address this 
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Fig. 23. — Sample overlap: U — BV and BV — Rci' colors for BzK galaxies (left) and z ~ 2 LBGs (right). Sources shown as black are 
those with Rc < 25.5 while those in gray are fainter than Rq = 25.5. Triangle symbols indicate non-detections in U. 
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Fig. 24. — Sample overlap: NUV — B and B — V colors for BX-, BM-, and BzK-selected galaxies. We have limited these sources to those 
with .Zphot ^ 1-5 since the NUV data are only able to see the Lyman continuum break at z ~ 1.5. Sources below V = 25.4 are shown as 
gray squares. Triangles are sources that are undetected in the NUV. 
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Fig. 25. — Accumulated SFR density against K- (top) and Rq- 
band (bottom) magnitudes from a census of star-forming galaxies 
at 1.5 < Zphot < 2.6. Shaded regions include uncertainties in the 
SFR due to the z p hot accuracy and cosmic variance estimates from 
Somervillc ct al. (2004). The census total is shown by the shaded 
gray region while those determined using the sBzK, BX/BM, and 
iVfTV-dropout methods are shown by the hatched red, green, and 
blue regions, respectively. A stellar mass or observed 1700 A SFR 
scale is provided for guidance. The stellar mass is from Daddi ct 
al. (2004) SED-fit relation. For the UV SFR, we adopt Kennicutt 
(1998) relation. (A color version of this figure is available in the 
online journal.) 

loss. 

7. MODELED PREDICTIONS FOR A MOCK CENSUS 

Motivation. To further understand the selection ef- 
fects of the census derived from color selection techniques 
(hereafter the "Census"), we generate Monte Carlo real- 
izations of our data and of the Census using mock galax- 
ies that represent the z « 1-3 population. This simula- 
tion intends to (1) examine the stellar population incom- 
pleteness effects on the multi-band x 2 -weighting scheme 



(discussed in Section 3.1), (2) create a census of mock 
LBGs, BX/BM, and BzK galaxies at z = 1.0-3.5, and 
(3) compare modeled predictions with observations. In 
summary, these analyses find both agreements and dis- 
agreements with observational results. The former firmly 
supports many of our main results while the discrepan- 
cies appear to have a minor effect. 

The ability for a galaxy to be included in the Cen- 
sus is dependent on a few factors. First, galaxies must 
be identified in the multi-band \ 2 image. We required 
X > 3.5 for at least five connecting pixels. Because the 
multi-band image spans a wide range in observed wave- 
lengths (4500 A-2 /xm), these x values are influenced by 
the shape of the SED, which can be parameterized by 
rcdshift, stellar mass, dust reddening, galaxy ages, and 
star formation histories. It is thus important to simulate 
a wide range of these galaxy properties (see below). Sec- 
ond, we place restrictions on the sample such that they 
are bright enough to be included in either the 5BK3cr or 
5B5a photo-z catalogs. Finally, we also apply the dif- 
ferent color selection techniques to cull LBGs, BX/BM, 
and BzK galaxies to generate a "Mock Census." These 
steps that we follow are identical to those conducted in 
Sections 3.1, 4.1, and 5. 

Technique/approach. We begin with a grid of spectral 
synthesis models that adopt a range of magnitudes, red- 
shifts, E(B — V) values, galaxy ages, and r values for 
an exponentially declining star formation history. These 
models allow us to generate the full SED for artificial 
galaxies, which are then convolved with the bandpasscs 
for all 15 bands to obtain apparent magnitudes. We then 
add noise based on known sensitivity at each bandpass to 
determine the probability of (1) being included into the 
multi-band x 2 image, (2) satisfying the 5BK3tr and/or 
5B5tr photo- z criteria, and (3) meeting the color selection 
criteria for the Census. 

Assumptions. In estimating the x values for mock 
galaxies, we assume that sources are unresolved with 
an FWHM of l'.'l. This process involves normalizing a 
Gaussian PSF by the S/N of the source in each band, and 
then taking the square-root of the sum of the squares for 
eight bands (BVR c i'z'K, IA598, and IA679). We then 
determine if the number of connecting 0'.'2 pixels that 
meets the minimum x = 3.5 criteria is at least five to be 
included. 

The grid of models spans redshifts from 1.0 to 3.2 (0.2 
increments), E(B — V) = 0.0-0.5 (0.1 increments), and 
ages between log(age/yr) = 7.0 dex and the age of the 
universe at a given redshift 11 (0.25 dex increments). We 
adopt exponentially declining star formation history with 
t = 0.1, 1.0, and 10.0 Gyr. The limits on these proper- 
ties are identical to those assumed in modeling the SEDs 
(see Section 4.2). We also assume a Chabrier (2003) IMF 
and solar metallicity for the spectral synthesis models. 
We also include redshift-dependent IGM absorption from 
neutral hydrogen for rest wavelengths blueward of 1216 
A following Madau (1995). With these parameters, we 
have a total of 2,268 SEDs. These SEDs are normalized 
in the z'-band between 20.0 and 26.0 mag (0.5 mag incre- 
ments), which yields 27,216 possible models. We chose to 
normalize at z' since it is the most sensitive band closest 

11 For 2 = 1 and 3, the limits are 5.75 and 2.11 Gyr, respectively. 
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Fig. 26. — Comparison between galaxies from color selections and the photo-z sample. The total numbers are shown on the left and the 
distributions are normalized to the photo-z sample (thin black solid lines) on the right. The sBzK, BX/BM, TVC/V-dropout, and [/-dropout 
samples are indicated by the dotted (red), short-dashed (green), dot-dashed (blue) and long-dashed lines (blue). The combination of these 
techniques, accounting for overlap between techniques, is shown by the thick gray solid lines. (A color version of this figure is available in 
the online journal.) 



to K. The use of the z'-band also provides a good com- 
promise between the rest-frame UV (sensitive to recent 
star formation) and optical (sensitive to stellar mass). 

We must adopt distributions for each modeled prop- 
erties (e.g., E(B — V), redshift, and magnitude). First, 
we assume that the observed z ~ 1-3 E(B — V) dis- 
tribution applies to all galaxies in the simulation. This 
simplifying assumption is needed for a proper normaliza- 
tion of the number of blue-to-red galaxies, since it will 
affect the overlap fraction between sBzK and BX/BM or 
LBGs and the fraction of the Census that is contributed 
by either method. The E{B — V) distribution shows a 
log-normal decline with higher dust reddening following 
log (TV) oc -Z.ZxE{B-V). Then for each E{B-V) bin, 
we use observational constraints on the redshift and mag- 
nitude distributions. The observed distributions show 
that the more reddened galaxies are found at z < 2 and 
have brighter z' magnitudes (see Figure 27). In total, 
we generate 34,755 mock galaxies and the simulation is 
repeated 100 times for 3,475,500 artificial galaxies. 

One caveat with these assumptions is that observed 
trends and distributions are used, which themselves are 
affected by incompleteness. For example, we could be 
underestimating the intrinsic number of extremely dusty 
and/or old galaxies. However, neither we nor previous 
investigators have sufficient information to attempt a fur- 
ther correction for this possible incompleteness. 

Predictions from simulations. With the Monte Carlo 
realization of our data, we find that the x 2 method begins 
to miss 5% (50%) of galaxies at V « 26.25 (V 27.25) 
and Rc ~ 25.75 (i?c ~ 26.5). This simulation shows 
that the % 2 method is highly complete down to mag- 
nitudes, V = 25.4 and Rc = 25.5, at which we select 
-/VCV-dropouts and BX/BM galaxies using optical col- 
ors. 

In Figure 28, we illustrate the fraction recovered by 



the x method as a function of dust reddening and the z'- 
band magnitude for the different r models. This figure 
illustrates high (>80%) completeness for the r = 1 and 
10 Gyr models regardless of brightness and E{B — V). As 
expected, the model with the shortest timescale for star 
formation will suffer more incompleteness since many of 
these galaxies are evolved and dust reddening further 
hampers their inclusion in the x 2 image. By coming all 
three r models, the incompleteness does not exceed 20%. 

Comparisons with observations. In addition to the 
above completeness estimates, our Monte Carlo simula- 
tion provides a few commensurable predictions. For ex- 
ample, in Figure 29 we illustrate the fraction of the Mock 
and Observed Census that is probed by each color selec- 
tion technique. This is similar to Figure 18 but with fixed 
redshift bins. There are several similarities that are ap- 
parent. The simulated ?7-dropout and BX distributions 
are fairly consistent with observations in the terms of the 
location and strength of the peak and the shape of the 
distributions. The sBzK distributions peak at z ~ 1.5 at 
80% (mock) and 70% (observed) and declines to 30%- 
40% at z > 3. Above z = 1.7, the z ~ 2 LBG distri- 
butions are similar peaking at z ~ 2 and extending out 
to z ~ 2.8. These similarities strongly suggest that the 
overlap determined between each population is a fairly 
reliable result, and it is indeed true that no technique 
yields 100% of the Census at any redshift. The greatest 
discrepancy, is at z < 1.5 where the mock simulation in- 
dicates higher completeness. The causes for such discrep- 
ancies is not known. However, these differences do not 
affect much of our main results which concern z > 1.5. 

Another prediction from the simulation is the fraction 
of the photo-z census that is recovered by the Census 
(see Figure 30). As we have shown, the observed census 
yields 90% of galaxies above z = 1.5, and this appears 
to hold out to z ~ 3. The Mock Census shows a similar 
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2.0 2.5 
Photometric Redshift 

Fig. 27. — Redshift (left) and .z'-band (right) distributions for galaxies with different measured E(B — V). These distributions are used 
as prior inputs for the Monte Carlo simulation. This figure illustrates that (observationally) most of the dusty galaxies are found at z < 2 
and at relatively brighter z' magnitudes. (A color version of this figure is available in the online journal.) 

behavior, but indicates nearly 100%. These differences 
are small since Poisson statistics is at the ~5% level. 

Finally, another set of positive comparisons is the V-, 
Rc~, and if-band luminosity functions for z ~ 2 LBGs. 
We find good to reasonable agreement over 5 magnitudes 
in the optical bands and 3 magnitudes in K . Other color 
selection techniques find agreements over a smaller range 
of magnitudes. The causes of poorer agreement are still 
unknown. 
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Fig. 28. — Recovered fraction of mock galaxies as a function of the 
assumed star formation history, dust reddening, and magnitude. 
Each panel plots E(B — V) as a function of z' . Results for r = 
0.1, 1, and 10 Gyr and all three combined are shown from top to 
bottom. The gray scale shows the completeness from 100% (white) 
to 0% (black). Faint magnitude bins with high reddening have large 
statistical fluctuations due to input small numbers. 



8. CONCLUSIONS 

We have conducted a large photometric survey of 
galaxies at z = 1-3 by synthesizing measurements from 
20 broad-, intermediate-, and narrow-band filters cover- 
ing observed wavelengths of 1500 A-4.5 (im. This sur- 
vey is unique due to its size (0.25 deg 2 ; 5.6 x 10 6 Mpc 3 
for z = 1~3), depth, and reliable photo-z's derived for 
w23,000 .fT-band selected sources («65,000 including op- 
tically selected sources). Compared to a spectroscopic 
sample, we find that our photo-z's from 20-bands are ac- 
curate to 1.2% x (1 + z) (la) for z < 1.8. The accuracy 
of our photo-z stems from coverage of both the Lyman 
continuum and Balmer/4000 A breaks at z ~ 2 with the 
multiple bands that we have. Most photo-z surveys use 
10 photometric bands or less, are limited to accurate red- 
shifts below z < 1 and z > 3, and suffer from a common 
photo-z problem of ambiguity between the Balmer/ 4000 
A break at z ~ 0.3 and Lyman break at z ~ 3. 

We have also examined the z p hot for «4500 narrow- 
band excess emitters at z « 0.25-1.5, and found excel- 
lent to adequate agreement between our derived z p hot, 
and the predictions from using simple multiple rest-frame 
optical colors (see Ly et al. 2007). This evidence further 
indicates that the photo-z that we have generated are 
reliable. 

In addition, the combination of deep optical and rela- 
tively deep NIR data provides a large and representative 
sample of z = 1-3 galaxies covering a wide range in stel- 
lar mass, dust content, and SFR. With ss21,000 («19,000 
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Fig. 29. — Comparison between the Mock Census (left) and the observed Census (right). The vertical axes show the fraction of the 
Census consisting of LBGs (black), BX (green), BM (blue), and sBzK (red) galaxies. These plots are a reproduction of Figure 18. Our 
Mock Census is able to reproduce the redshift peak of each galaxy population and the shape of the redshift distributions above z ~ 1.5. 
These comparisons are further discussed in Section 7. (A color version of this figure is available in the online journal.) 

This problem is largely attributed to the Balmcr/4000 
A break occurring between the U and B filters, rather 
than a Lyman break. We also find that the NUV-dropout 
technique has 17.4% contamination from "interlopers" 
at z = 1.0-1.5, which eliminates a previous concern that 
contamination fraction estimates in Ly et al. (2009) were 
too optimistic. Also, the technique does not suffer from 
z < 0.5 interlopers, since the Balmer break is rcdward of 
the NUV. 

Selection bias and completeness of techniques. We con- 
firmed that each color selection is biased towards some 
subset of galaxies not fully representative of the entire 
galaxy population. By modeling the SED of thousands of 
photometrically-selected galaxies, we find that the BzK 
method targets more massive and dusty galaxies, while 
the BX, BM, and [/-dropout methods identify galaxies 
of lower stellar mass with low reddening. Compared to 
a combined census of color-selected galaxies (BX/BM, 
LBG, and BzK), we find that any UV or IR technique 
finds at most 60%-80% at the peak of their redshift sensi- 
tivity. However, these techniques arc optimized for a cer- 
tain redshift range, so relative toaz = 1-3 color-selected 
census, the BX/BM, ./VCV-dropout and sBzK techniques 
obtain 52%-65%. This strongly indicates that combin- 
ing multiple techniques is favored for a complete census 
of high-redshift galaxies and their SFRs. We estimate 
that the BX/BM method is capable of identifying (by 
number) 70% of z p hot = 1 — 3 census galaxies with SFR 
between a few and 30 M Q yr _1 , and this represents 66% 
of the Zphot = 1~3 census. 

Sample overlap between galaxy populations. Among 
the star-forming BzK-sclected galaxies, the overlap frac- 
tion with the BX/BM/ [/-dropout increases at fainter K 
luminosities, consistent with R05. However, with a larger 
sample of massive sBzK galaxies, we determined that the 
overlap fraction is ~5 times smaller than what was re- 
ported in R05 at bright K magnitudes. We argue that 
this is due to the smaller FoV and sample size that they 
used. Our study also shows that 77.2% of z > 1.5 sBzK 
galaxies have a strong Lyman break to meet the NUV- 
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Fig. 30. — Comparison between the Mock Census (solid line) and 
the observed Census (dashed line). The vertical axis shows the 
fraction of the photo-z sample that consists of galaxies that satisfy 
either one of the criteria for being BX, BM, sBzK and LBGs at 
z ~ 2 and z ~ 3. This is a reproduction of Figure 26 (right). The 
Mock Census yields a higher fraction than observed; however, it 
is consistent to within 2cr of Poisson statistics. These comparisons 
are further discussed in Section 7. 

with Zphot) galaxies, we are able to study statistically the 
overlap of different galaxy populations (BX/BM, LBG, 
and BzK) obtained from each technique. We also in- 
vestigate the completeness of these techniques against a 
sample derived purely from photo-z. The main results 
for this paper are: 

Redshift distribution. We determined (with accurate 
photo-z) that 74.7% of BXs, 84.5% of BMs, 78.83% of U- 
dropouts, 93.9% of sBzKs, and 91.9% of A/W-dropouts 
have Zphot = 1-3.5. We find that the BX, BM, and 
{/-dropout methods suffer from z p hot < 0.5 contamina- 
tion at the 24.4%, 12.8%, and 13.4% levels, respectively. 
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dropout criteria. Unlike the BX/BM method, this over- 
lap fraction is independent of K mag, indicating that 
the Lyman break technique is relatively more sensitive 
to massive galaxies. Likewise, we find that 79.7% of 
z > 1.5 BX and 87.0% of z > 1.5 BM galaxies are z ~ 2 
LBGs. These results indicate that the 912 A break is 
ubiquitous in z ~ 2 galaxies. Thus, the Af[/U-dropout 
selection differs from the BX/BM selection by including 
more massive and reddened galaxies. 

SFR density from a census of UV- and NIR-selected 
galaxies. We determined that the comoving SFR density 
at 2 phot = 1.5-2.6 is 0.18±0.03 M Q yr" 1 Mpc~ 3 from 
a joint census of BzK galaxies, BX/BM galaxies, and 
LBGs. This census was limited to galaxies with K < 24 
(corresponds to « 10 10 M®). We find that the UV selec- 
tion methods (Lyman break and BX/BM) obtain 81% of 
the SFR density from galaxies with E(B — V) < 0.25 
mag. These low reddening galaxies represent about 
three-quarters of our census by number. However, there 
is a prominent population of galaxies with high dust ex- 
tinction (E(B — V) > 0.25 mag). For example, we ex- 
amined the contribution to the SFR density relative to 
the dust properties of our galaxies, and find that galaxies 
below and above E[B — V) = 0.25 mag contribute 35% 
and 65% to the SFR density, respectively. 

Completeness of color selection techniques. We com- 
pare the census derived from color selection techniques 
against a photo-z sample, and find that color selection 
is (1) less efficient at z < 1.5, but (2) identifies «90% 
of galaxies at z p hot = 1.5-2.5. We also find that the 
sBzK method is 85% complete when compared against 
a sample of if -band detected sources with z p hot > 1.5. 
This evidence, combined with the high fraction of over- 
lap between BX/BM and sBzK galaxies for K ~ 23-24 
AB mag, and the likelihood that faint (blue) UV-selected 



galaxies will be detected with deeper K-h&nd imaging, 
strongly suggest that future z « 1-3 studies should focus 
on a very deep A-band imaging (e.g., the Ultra- VISTA 
survey). But until such a deep (A ~ 25-26 AB) survey 
with multi-wavelength data is available over large areas, 
the UV selection techniques (cither the U n GlZ or NUV- 
dropout method) are currently the only efficient way to 
identify less massive z ~ 2 galaxies. 
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APPENDIX 

A. UNUSUAL SOURCES FROM COLOR SELECTION 

About half of the Shallow and Faint samples (restricted to the deep A-band area and excluding galaxies with 
Zphot < 0.5) consists of sources that are uniquely selected by only one color selection. These sources arc of great interest, 
since they indicate that a particular technique culls a galaxy sample that is unavailable by any other techniques. We 
summarize the uniqueness below. 

Passive BzK galaxies. We find that the entire pBzK sample does not overlap with any of the other techniques. These 
galaxies at z sw 1-2 have the oldest ages, log (t age /yr) ss 9.0, the highest stellar masses, log (M/M Q ) sa 11.0, and span 
E(B — V) = 0.0-0.3. Their ages imply that they formed Az = 0.6 earlier (median formation redshift is z = 2.25), and 
could be the descendants of luminous infrared galaxies or massive LBGs with SFRs > 100 M Q yr" 1 . These galaxies 
are important for a stellar mass census survey, but contribute little to the SFR density since their SFRs are on average 
0.02 M Q yr- 1 . 

Star-forming BzK galaxies. We find that purely sBzKs, compared to UV-selected samples, have similar SFRs (~ 5 
Mq yr -1 ), but are more massive (~ 10 9,5 -10 11 M Q ), cover a broader range in dust extinction (E[B — V] « 0.05-0.4), 
and are 0.5 dex older. 

BX, BM, and U -dropout galaxies. The UV selection techniques of Steidel et al. identify a unique set of galaxies 
with relatively younger stellar ages (typically 2 x 10 8 yr), lower stellar mass (~ 10 9 Mq), and E(B — V) rts 0.1. We 
find that the SFRs are systematically lower by 0.2-0.3 dex for BM and BX galaxies when compared to [/-dropouts, 
although this is likely a manifestation of redshift dependence of the sample. The relatively poor [/-dropout overlap 
with NIR techniques can be attributed to the redshift of the sample, as well as, the lack of deeper NIR data. It would 
be prudent to identify a sample of higher redshift BzK galaxies using a set of different filters (e.g., R, J, and L; Daddi 
et al. 2004) to investigate the overlap of LBGs at higher redshift with 3.6 /zm-selected galaxies. 

Lyman break galaxies at z ~ 2. The method that we developed to identify LBGs at z sa 1.5-2.5 has significant 
overlap («80%) with the BX/BM and/or BzK techniques. It appears that the remaining 20% consists of galaxies 
with stellar masses of « 10 10±1 M Q , cover the same stellar ages as other techniques, have an average SFR of ^20 M 
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yr^ 1 , and span E(B — V) = 0.1-0.4. The main reason for why these galaxies were uniquely identified is their redshift: 
about two-thirds of them are below z — 1.5 where the sBzK and BM techniques are less sensitive. We note that the 
higher SFRs for A^t/V-dropouts are likely due to the limit adopted for the X^-band (our limit is comparable brighter 
than Rc = 25.5 for BX/BM). The fainter (V = 26.0) limit yields an average SFR that is systematically lower by 0.2 
dex, due to the higher contribution from lower-SFR LBGs. 
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